
> #########################################################################
> # 	SCRIPT FOR "RECONSIDERING THE NEIGHBORHOOD EFFECT"			#
> #	Authors: Martin Bisgaard, Peter T. Dinesen & Kim S�nderskov 	#
> #	Last edited: 10-11-2015
> 
> 
> rm(list=ls())
> require(foreign)
> #load data
> dat<-read.dta("mydirectory")
> names(dat)
  [1] "pnr"              "idno"             "essround"        
  [4] "panel"            "idate"            "idate2"          
  [7] "intwqtr1"         "intwqtr2"         "dayslived"       
 [10] "dayslivedafter"   "waves"            "move"            
 [13] "stayed"           "moved3after"      "knr"             
 [16] "kom"              "municip"          "stfeco"          
 [19] "stfgov"           "feelhhincome"     "lifesat"         
 [22] "lrscale"          "prtvtdk2002"      "prtvtdk2004"     
 [25] "prtvtdk2006"      "prtvtdk2008"      "prtvtdk2010"     
 [28] "prtvtdk2011"      "polintr"          "polcmpl"         
 [31] "poldcs"           "polmedia"         "tvtot"           
 [34] "tvpol"            "rdtot"            "rdpol"           
 [37] "nwsptot"          "nwsppol"          "mainactyear"     
 [40] "mainactweek"      "unemplo2"         "unemplo3"        
 [43] "unemplo4"         "unemplo5"         "pincome2000"     
 [46] "pincome2002"      "pincome2004"      "pincome2006"     
 [49] "pincome2008"      "pincome2010"      "pincome2011"     
 [52] "shareunem"        "shareimm"         "noimbro"         
 [55] "ST"               "soctru1"          "itru"            
 [58] "maleR"            "maleS"            "maleRS"          
 [61] "ageR"             "ageRS"            "partunemplo4"    
 [64] "oneunemplo4"      "flag"             "cunemploy"       
 [67] "cunmminus1y"      "munemploy"        "munemploym1y"    
 [70] "edusur"           "edureg"           "eduresu"         
 [73] "immdk"            "nww1g"            "nww2g"           
 [76] "wdknw1g"          "wdknw2g"          "muslim"          
 [79] "muslim2"          "immdk80"          "immdk130"        
 [82] "immdk180"         "immdk250"         "immdk350"        
 [85] "immdk500"         "immdk750"         "immdk1000"       
 [88] "immdk1250"        "immdk1500"        "immdk1750"       
 [91] "immdk2000"        "immdk2250"        "immdk2500"       
 [94] "nww1g80"          "nww2g80"          "nww1g130"        
 [97] "nww2g130"         "nww1g180"         "nww2g180"        
[100] "nww1g250"         "nww2g250"         "nww1g350"        
[103] "nww2g350"         "nww1g500"         "nww2g500"        
[106] "nww1g750"         "nww2g750"         "nww1g1000"       
[109] "nww2g1000"        "nww1g1250"        "nww2g1250"       
[112] "nww1g1500"        "nww2g1500"        "nww1g1750"       
[115] "nww2g1750"        "nww1g2000"        "nww2g2000"       
[118] "nww1g2250"        "nww2g2250"        "nww1g2500"       
[121] "nww2g2500"        "wdknw1g80"        "wdknw2g80"       
[124] "wdknw1g130"       "wdknw2g130"       "wdknw1g180"      
[127] "wdknw2g180"       "wdknw1g250"       "wdknw2g250"      
[130] "wdknw1g350"       "wdknw2g350"       "wdknw1g500"      
[133] "wdknw2g500"       "wdknw1g750"       "wdknw2g750"      
[136] "wdknw1g1000"      "wdknw2g1000"      "wdknw1g1250"     
[139] "wdknw2g1250"      "wdknw1g1500"      "wdknw2g1500"     
[142] "wdknw1g1750"      "wdknw2g1750"      "wdknw1g2000"     
[145] "wdknw2g2000"      "wdknw1g2250"      "wdknw2g2250"     
[148] "wdknw1g2500"      "wdknw2g2500"      "muslim80"        
[151] "muslim130"        "muslim180"        "muslim250"       
[154] "muslim350"        "muslim500"        "muslim750"       
[157] "muslim1000"       "muslim1250"       "muslim1500"      
[160] "muslim1750"       "muslim2000"       "muslim2250"      
[163] "muslim2500"       "popdens80"        "popdens130"      
[166] "popdens180"       "popdens250"       "popdens350"      
[169] "popdens500"       "popdens750"       "popdens1000"     
[172] "popdens1250"      "popdens1500"      "popdens1750"     
[175] "popdens2000"      "popdens2250"      "popdens2500"     
[178] "numobs1_80"       "numobs1_130"      "numobs1_180"     
[181] "numobs1_250"      "numobs1_350"      "numobs1_500"     
[184] "numobs1_750"      "numobs1_1000"     "numobs1_1250"    
[187] "numobs1_1500"     "numobs1_1750"     "numobs1_2000"    
[190] "numobs1_2250"     "numobs1_2500"     "ethdiv80"        
[193] "ethdiv130"        "ethdiv180"        "ethdiv250"       
[196] "ethdiv350"        "ethdiv500"        "ethdiv750"       
[199] "ethdiv1000"       "ethdiv1250"       "ethdiv1500"      
[202] "ethdiv1750"       "ethdiv2000"       "ethdiv2250"      
[205] "ethdiv2500"       "aincome2002_80"   "aincome2000_80"  
[208] "aincome2002_130"  "aincome2000_130"  "aincome2002_180" 
[211] "aincome2000_180"  "aincome2002_250"  "aincome2000_250" 
[214] "aincome2002_350"  "aincome2000_350"  "aincome2002_500" 
[217] "aincome2000_500"  "aincome2002_750"  "aincome2000_750" 
[220] "aincome2002_1000" "aincome2000_1000" "aincome2002_1250"
[223] "aincome2000_1250" "aincome2002_1500" "aincome2000_1500"
[226] "aincome2002_1750" "aincome2000_1750" "aincome2002_2000"
[229] "aincome2000_2000" "aincome2002_2250" "aincome2000_2250"
[232] "aincome2002_2500" "aincome2000_2500" "uen80"           
[235] "uen130"           "uen180"           "uen250"          
[238] "uen350"           "uen500"           "uen750"          
[241] "uen1000"          "uen1250"          "uen1500"         
[244] "uen1750"          "uen2000"          "uen2250"         
[247] "uen2500"          "ueb80"            "ueb130"          
[250] "ueb180"           "ueb250"           "ueb350"          
[253] "ueb500"           "ueb750"           "ueb1000"         
[256] "ueb1250"          "ueb1500"          "ueb1750"         
[259] "ueb2000"          "ueb2250"          "ueb2500"         
[262] "numobs2_80"       "numobs2_130"      "numobs2_180"     
[265] "numobs2_250"      "numobs2_350"      "numobs2_500"     
[268] "numobs2_750"      "numobs2_1000"     "numobs2_1250"    
[271] "numobs2_1500"     "numobs2_1750"     "numobs2_2000"    
[274] "numobs2_2250"     "numobs2_2500"     "numobs3_80"      
[277] "numobs3_130"      "numobs3_180"      "numobs3_250"     
[280] "numobs3_350"      "numobs3_500"      "numobs3_750"     
[283] "numobs3_1000"     "numobs3_1250"     "numobs3_1500"    
[286] "numobs3_1750"     "numobs3_2000"     "numobs3_2250"    
[289] "numobs3_2500"     "aincome2004_80"   "aincome2004_130" 
[292] "aincome2004_180"  "aincome2004_250"  "aincome2004_350" 
[295] "aincome2004_500"  "aincome2004_750"  "aincome2004_1000"
[298] "aincome2004_1250" "aincome2004_1500" "aincome2004_1750"
[301] "aincome2004_2000" "aincome2004_2250" "aincome2004_2500"
[304] "aincome2006_80"   "aincome2006_130"  "aincome2006_180" 
[307] "aincome2006_250"  "aincome2006_350"  "aincome2006_500" 
[310] "aincome2006_750"  "aincome2006_1000" "aincome2006_1250"
[313] "aincome2006_1500" "aincome2006_1750" "aincome2006_2000"
[316] "aincome2006_2250" "aincome2006_2500" "aincome2008_80"  
[319] "aincome2008_130"  "aincome2008_180"  "aincome2008_250" 
[322] "aincome2008_350"  "aincome2008_500"  "aincome2008_750" 
[325] "aincome2008_1000" "aincome2008_1250" "aincome2008_1500"
[328] "aincome2008_1750" "aincome2008_2000" "aincome2008_2250"
[331] "aincome2008_2500" "aincome2010_80"   "aincome2010_130" 
[334] "aincome2010_180"  "aincome2010_250"  "aincome2010_350" 
[337] "aincome2010_500"  "aincome2010_750"  "aincome2010_1000"
[340] "aincome2010_1250" "aincome2010_1500" "aincome2010_1750"
[343] "aincome2010_2000" "aincome2010_2250" "aincome2010_2500"
[346] "aincome2011_80"   "aincome2011_130"  "aincome2011_180" 
[349] "aincome2011_250"  "aincome2011_350"  "aincome2011_500" 
[352] "aincome2011_750"  "aincome2011_1000" "aincome2011_1250"
[355] "aincome2011_1500" "aincome2011_1750" "aincome2011_2000"
[358] "aincome2011_2250" "aincome2011_2500" "aincome2yb80"    
[361] "aincome2yb130"    "aincome2yb180"    "aincome2yb250"   
[364] "aincome2yb350"    "aincome2yb500"    "aincome2yb750"   
[367] "aincome2yb1000"   "aincome2yb1250"   "aincome2yb1500"  
[370] "aincome2yb1750"   "aincome2yb2000"   "aincome2yb2250"  
[373] "aincome2yb2500"   "uen2yb80"         "uen2yb130"       
[376] "uen2yb180"        "uen2yb250"        "uen2yb350"       
[379] "uen2yb500"        "uen2yb750"        "uen2yb1000"      
[382] "uen2yb1250"       "uen2yb1500"       "uen2yb1750"      
[385] "uen2yb2000"       "uen2yb2250"       "uen2yb2500"      
[388] "ueb2yb80"         "ueb2yb130"        "ueb2yb180"       
[391] "ueb2yb250"        "ueb2yb350"        "ueb2yb500"       
[394] "ueb2yb750"        "ueb2yb1000"       "ueb2yb1250"      
[397] "ueb2yb1500"       "ueb2yb1750"       "ueb2yb2000"      
[400] "ueb2yb2250"       "ueb2yb2500"       "aincome2ya80"    
[403] "aincome2ya130"    "aincome2ya180"    "aincome2ya250"   
[406] "aincome2ya350"    "aincome2ya500"    "aincome2ya750"   
[409] "aincome2ya1000"   "aincome2ya1250"   "aincome2ya1500"  
[412] "aincome2ya1750"   "aincome2ya2000"   "aincome2ya2250"  
[415] "aincome2ya2500"   "uen2ya80"         "uen2ya130"       
[418] "uen2ya180"        "uen2ya250"        "uen2ya350"       
[421] "uen2ya500"        "uen2ya750"        "uen2ya1000"      
[424] "uen2ya1250"       "uen2ya1500"       "uen2ya1750"      
[427] "uen2ya2000"       "uen2ya2250"       "uen2ya2500"      
[430] "ueb2ya80"         "ueb2ya130"        "ueb2ya180"       
[433] "ueb2ya250"        "ueb2ya350"        "ueb2ya500"       
[436] "ueb2ya750"        "ueb2ya1000"       "ueb2ya1250"      
[439] "ueb2ya1500"       "ueb2ya1750"       "ueb2ya2000"      
[442] "ueb2ya2250"       "ueb2ya2500"       "singchilhh80"    
[445] "singlehh80"       "singchilhh130"    "singlehh130"     
[448] "singchilhh180"    "singlehh180"      "singchilhh250"   
[451] "singlehh250"      "singchilhh350"    "singlehh350"     
[454] "singchilhh500"    "singlehh500"      "singchilhh750"   
[457] "singlehh750"      "singchilhh1000"   "singlehh1000"    
[460] "singchilhh1250"   "singlehh1250"     "singchilhh1500"  
[463] "singlehh1500"     "singchilhh1750"   "singlehh1750"    
[466] "singchilhh2000"   "singlehh2000"     "singchilhh2250"  
[469] "singlehh2250"     "singchilhh2500"   "singlehh2500"    
[472] "conedu80"         "conedu130"        "conedu180"       
[475] "conedu250"        "conedu350"        "conedu500"       
[478] "conedu750"        "conedu1000"       "conedu1250"      
[481] "conedu1500"       "conedu1750"       "conedu2000"      
[484] "conedu2250"       "conedu2500"       "congymreg80"     
[487] "congymreg130"     "congymreg180"     "congymreg250"    
[490] "congymreg350"     "congymreg500"     "congymreg750"    
[493] "congymreg1000"    "congymreg1250"    "congymreg1500"   
[496] "congymreg1750"    "congymreg2000"    "congymreg2250"   
[499] "congymreg2500"    "Dctlvlumemp"      "Dmulvlunemp"     
[502] "Duen80"           "Duen130"          "Duen180"         
[505] "Duen250"          "Duen350"          "Duen500"         
[508] "Duen750"          "Duen1000"         "Duen1250"        
[511] "Duen1500"         "Duen1750"         "Duen2000"        
[514] "Duen2250"         "Duen2500"         "Dueb80"          
[517] "Dueb130"          "Dueb180"          "Dueb250"         
[520] "Dueb350"          "Dueb500"          "Dueb750"         
[523] "Dueb1000"         "Dueb1250"         "Dueb1500"        
[526] "Dueb1750"         "Dueb2000"         "Dueb2250"        
[529] "Dueb2500"         "Dinco80"          "Dinco130"        
[532] "Dinco180"         "Dinco250"         "Dinco350"        
[535] "Dinco500"         "Dinco750"         "Dinco1000"       
[538] "Dinco1250"        "Dinco1500"        "Dinco1750"       
[541] "Dinco2000"        "Dinco2250"        "Dinco2500"       
[544] "dup"             

> dim(dat)
[1] 7684  544

> ####################################
> # RESCALING AND RECODING VARIABLES #
> 
> require(car)

> #main DV: economic perceptions
> dat$stfeco<-dat$stfeco/10 #rescaling to 0-1

> table(dat$stfeco) 

   0  0.1  0.2  0.3  0.4  0.5  0.6  0.7  0.8  0.9    1 
  93   64  215  385  496  873  863 1384 1775  839  483 

> #self-reported exposure to TV, Radio and newspapers
> dat$tvtot<-recode(dat$tvtot,'88=NA;99=NA');dat$tvtot<-dat$tvtot/7

> table(dat$tvtot)

                0 0.142857142857143 0.285714285714286 0.428571428571429 
              116               377              1135              1210 
0.571428571428571 0.714285714285714 0.857142857142857                 1 
             1347              1138               971              1382 

> dat$rdtot<-recode(dat$rdtot,'88=NA;99=NA');dat$rdtot<-dat$rdtot/7

> table(dat$rdtot)

                0 0.142857142857143 0.285714285714286 0.428571428571429 
             1095              1506              1239               628 
0.571428571428571 0.714285714285714 0.857142857142857                 1 
              415               329               268              2191 

> dat$nwsptot<-recode(dat$nwsptot,'88=NA;99=NA');dat$nwsptot<-dat$nwsptot/7

> table(dat$nwsptot)

                0 0.142857142857143 0.285714285714286 0.428571428571429 
             1802              2682              2203               622 
0.571428571428571 0.714285714285714 0.857142857142857                 1 
              208               101                36                25 

> #satisfaction with government (used for supplemental analyses)
> dat$stfgov<-dat$stfgov/10

> table(dat$stfgov)

   0  0.1  0.2  0.3  0.4  0.5  0.6  0.7  0.8  0.9    1 
 247  232  434  728  750 1101  942 1258 1209  407  174 

> #vote intention (used for supplemental analyses)
> dat$pvote<-NA

> dat$pvote[dat$essround == 1]<-dat$prtvtdk2002[dat$essround == 1]

> dat$pvote[dat$essround == 2]<-dat$prtvtdk2004[dat$essround == 2]

> dat$pvote[dat$essround == 3]<-dat$prtvtdk2006[dat$essround == 3]

> dat$pvote[dat$essround == 4]<-dat$prtvtdk2008[dat$essround == 4]

> dat$pvote[dat$essround == 5]<-recode(dat$prtvtdk2010[dat$essround == 5],'4=5;5=6;6=7;7=8;8=12;9=10;10=11')

> dat$pvote<-recode(dat$pvote,'66=13;77=13;88=13;99=13')

> table(dat$pvote)

   1    2    3    4    5    6    7    8    9   10   11   13 
1407  342  423  235  432  293  436 1288   41   91   17 1103 

> #recoding to vote intention to binary: incumbent parties/opposition parties
> dat$gov<-NA

> dat$gov[dat$pvote%in%c(3,8)]<-1

> dat$gov[dat$pvote%in%c(1:2,4:7,9:13)]<-0

> table(dat$gov)

   0    1 
4397 1711 

> # rescaling additional measures
> dat$munemploy<-dat$munemploy/100 #municipality unemployment

> dat$pincome2000_1k<-dat$pincome2000/100000 

> dat$dayslived_y<-dat$dayslived/365 #scaling to years lived at location

> dat$lrscale<-dat$lrscale/10 #left-right orientation

> ################################
> #  Code producing Figure 1	#
> 
> #Results plotted in Panel A
> #computing differences from 80m context
> 
> ueb.names<-paste("ueb",c(80,130,180,250,350,500,750,1000,1250,1500,1750,2000,2250,2500),sep="")

> diff<-dat[dat$popdens80>9,ueb.names]#matrix with all contextual measures

> diff.out<-matrix(NA,nr=nrow(dat[dat$essround!=9&dat$popdens80>9,]),nc=length(ueb.names)-1)

> for(i in 1:ncol(diff.out)) diff.out[,i]<-diff[,i+1]-diff[,1]

> #Results plotted in Panel B
> 
> 
> #inter-correlation between 80m and all other contexts (solid line in figure)
> 
> coef.out<-matrix(NA,nr=14,nc=2)#matrix to hold output

> for (i in 2:14){
+ form<-substitute(ueb80 ~ sub,list(sub=as.name(ueb.names[i]))) #model formula
+ fit1<-lm(form,data=dat[dat$popdens80>9,]) #model fit (restricted to respondents with more than 9 ppl in their 80m radius)
+ coef.out[i,1]<-coef(fit1)[ueb.names[i]]#extract coefficient
+ 	SSreg<-sum(fit1$residuals^2)/fit1$df
+ 	SStot<-var(fit1$model[,1])
+ coef.out[i,2]<-1-SSreg/SStot #extract adj. rsquared
+ }

> coef.out
           [,1]      [,2]
 [1,]        NA        NA
 [2,] 1.0106518 0.7577906
 [3,] 1.0050771 0.6393275
 [4,] 0.9934838 0.5281288
 [5,] 1.0001554 0.4388750
 [6,] 0.9726004 0.3416162
 [7,] 0.9540355 0.2619767
 [8,] 0.9539830 0.2194567
 [9,] 0.9437204 0.1890961
[10,] 0.9349725 0.1684184
[11,] 0.9333522 0.1538519
[12,] 0.9408210 0.1427553
[13,] 0.9464604 0.1314335
[14,] 0.9437753 0.1210485

> coef.out[1,]<-1 #setting an artificial '1' for 80m on 80m regression

> #inter-correlation between municipality unemployment and all other contexts (dashed line in figure)
> coef.out1<-matrix(NA,nr=14,nc=2)

> for (i in 1:14){
+ form<-substitute(munemploy ~ sub,list(sub=as.name(ueb.names[i])))
+ fit1<-lm(form,data=dat[dat$essround!=9&dat$popdens80>9,])
+ 
+ coef.out1[i,1]<-coef(fit1)[ueb.names[i]] #coef
+ 	SSreg<-sum(fit1$residuals^2)/fit1$df
+ 	SStot<-var(fit1$model[,1])
+ coef.out1[i,2]<-1-SSreg/SStot #adj. rsquared
+ }

> coef.out1
            [,1]       [,2]
 [1,] 0.04101298 0.03989643
 [2,] 0.05444138 0.05237261
 [3,] 0.06398297 0.06179088
 [4,] 0.07431865 0.07066057
 [5,] 0.08733793 0.08014382
 [6,] 0.10521363 0.09591971
 [7,] 0.13518309 0.12598092
 [8,] 0.15820747 0.14463387
 [9,] 0.17902932 0.16308911
[10,] 0.19775046 0.18053224
[11,] 0.21537369 0.19615716
[12,] 0.23430091 0.21179020
[13,] 0.25438195 0.22701205
[14,] 0.27158465 0.23961556

> #	PLOTTING FIGURE 1	#
> 
> #calculating the proportions of distribution below/over 5% (depicted in Panel A, Figure 1)
> below=prop.table(table(diff.out[,13]<c(-.05)))[2]

> above=prop.table(table(diff.out[,13]>c(.05)))[2]

> remain=1-(below+above)

> below;above;remain
     TRUE 
0.1935435 
     TRUE 
0.3668342 
     TRUE 
0.4396224 

> below+above+remain
TRUE 
   1 

> #pdf("bivar.pdf",width=7.75,height=4)
> setEPS()

> postscript("bivar.eps",width=7.75,height=4)

> xax<-c(80,130,180,250,350,500,750,1000,1250,1500,1750,2000,2250,2500)#spacing for x-axis (Panel B)

> par(mfrow=c(1,2))

> dens=density(diff.out[,13],na.rm=T,bw=0.005)

> plot(dens,xlim=c(-.3,.2),main="",xlab="",yaxt='n',ylab="Kernel density")

> axis(2,las=2)

> mtext(1,text="Deviation in unemployment rates",line=2.5)

> mtext(1,text="(2,500m - 80m)",line=3.5)

> mtext(3,text="A",font=2,cex=2,line=2)

> mtext(3,text="How much does aggregate unemployment (2,500m)",line=1,cex=.8)

> mtext(3,text="deviate from residential unemployment (80m)?",line=.25,cex=.8)

> mtext(2,text="(bandwidth = 0.005)",cex=.8,line=2)

> x1<-min(which(dens$x>=.05))

> x2<-max(which(dens$x<c(-.05)))

> with(dens,polygon(x=c(x[c(x1,x1:x2,x2)]),y=c(0,y[x1:x2],0),col="grey80"))

> text(y=0.35,x=-.1,label="19 %")

> text(y=0.35,x=0,label="44 %")

> text(y=0.35,x=.1,label="37 %")

> plot(x=xax,coef.out[,2],ylim=c(0,1),yaxt="n",xaxt="n",pch=1,col="white",
+ 	xlab="",ylab="R-Squared")

> mtext(3,text="B",font=2,cex=2,line=2)

> mtext(3,text="The proportion of variance explained",line=1,cex=.8)

> mtext(3,text="by different contextual measures",line=.25,cex=.8)

> lines(xax,coef.out[,2],lwd=1)

> lines(xax,coef.out1[,2],lwd=1,lty=2)

> axis(2,seq(0,1,.2),las=2,cex.axis=1)

> pos<-c(0,250,500,750,1000,1250,1500,1750,2000,2250,2500)

> axis(1,at=pos,labels=NA)

> text(pos[-1],par()$usr[3]-.06*(par()$usr[4]-par()$usr[3]),labels=c("250m","500m","750m","1000m","1250m","1500m","1750m","2000m","2250m","2500m"),
+ 	adj=1,srt=45,xpd=T,cex=.8)

> legend("topright",title="Dependent variable:",bty="n",legend=c("80m unemployment","Municipality unemployment"),lty=c(1,2),cex=.8,xpd=TRUE,lwd=1,pch=c(21,NA,1,NA),pt.bg="black")

> points(xax,coef.out[,2],pch=21,cex=1,bg="black")

> points(xax,coef.out1[,2],pch=1,cex=1,bg="black")

> mtext(1,text="Level of aggregation of contextual",line=2.5)

> mtext(1,text="predictor (radius in meters)",line=3.5)

> dev.off()
null device 
          1 

> ##########################
> # Results for Table 1	#
> require(memisc)
Indl�ser kr�vet pakke: memisc
Indl�ser kr�vet pakke: lattice
Indl�ser kr�vet pakke: MASS

Vedh�fter pakke: �memisc�

Det f�lgende objekt er maskeret fra �package:car�:

    recode

De f�lgende objekter er maskerede fra �package:stats�:

    contr.sum, contr.treatment, contrasts

De f�lgende objekter er maskerede fra �package:base�:

    as.array, trimws


> fit1<-lm(stfeco~factor(municip)+factor(essround)+popdens80+unemplo4+ageR+eduresu+maleR+pincome2000_1k+dayslived_y+lrscale+nwsptot+polintr+ueb80
+ 	,data=dat[dat$popdens80>9,])

> fit2<-lm(stfeco~factor(municip)+factor(essround)+popdens80+unemplo4+ageR+eduresu+maleR+pincome2000_1k+dayslived_y+lrscale+nwsptot+polintr+ueb80+ueb2500
+ 	,data=dat[dat$popdens80>9,])

> fit3<-lm(stfeco~factor(municip)+factor(essround)+popdens80+unemplo4+ageR+eduresu+maleR+pincome2000_1k+dayslived_y+lrscale+nwsptot+polintr+ueb80+ueb2500+oneunemplo4
+ 	,data=dat[dat$popdens80>9,])

> tab<-mtable(fit1,fit2,fit3)

> tab<-relabel(tab,
+ 	"ueb80"="Local Unemployment Rate (80m)",
+ 	ueb2500="Local Unemployment Rate (2500m)",
+ 	ueb2500="Local Unemployment Rate (2500m)",
+ 	oneunemplo4="Household Unemployment",
+ 	lrscale="Right--wing",
+ 	polintr="Political interest",
+ 	nwsptot="Newspaper exposure",
+ 	dayslived_y="Years lived at location",
+ 	pincome2000_1k="Income (in 100k DKR)",
+ 	ageR="Age",
+ 	unemplo4="Unemployed",
+ 	popdens80="Population (80m)",
+ 	eduresu="Years of education",
+ 	"maleR: Male/Female"="Male"
+ 	)

> tab

Calls:
fit1: lm(formula = stfeco ~ factor(municip) + factor(essround) + popdens80 + 
    unemplo4 + ageR + eduresu + maleR + pincome2000_1k + dayslived_y + 
    lrscale + nwsptot + polintr + ueb80, data = dat[dat$popdens80 > 
    9, ])
fit2: lm(formula = stfeco ~ factor(municip) + factor(essround) + popdens80 + 
    unemplo4 + ageR + eduresu + maleR + pincome2000_1k + dayslived_y + 
    lrscale + nwsptot + polintr + ueb80 + ueb2500, data = dat[dat$popdens80 > 
    9, ])
fit3: lm(formula = stfeco ~ factor(municip) + factor(essround) + popdens80 + 
    unemplo4 + ageR + eduresu + maleR + pincome2000_1k + dayslived_y + 
    lrscale + nwsptot + polintr + ueb80 + ueb2500 + oneunemplo4, 
    data = dat[dat$popdens80 > 9, ])

==============================================================
                                   fit1      fit2      fit3   
--------------------------------------------------------------
(Intercept)                       0.602***  0.598***  0.601***
                                 (0.023)   (0.028)   (0.028)  
factor(municip): 147/101          0.002     0.002     0.002   
                                 (0.021)   (0.021)   (0.021)  
factor(municip): 151/101         -0.076*   -0.076*   -0.074*  
                                 (0.037)   (0.037)   (0.037)  
factor(municip): 153/101          0.001     0.002     0.003   
                                 (0.039)   (0.039)   (0.039)  
factor(municip): 155/101         -0.018    -0.017    -0.018   
                                 (0.044)   (0.044)   (0.044)  
factor(municip): 157/101         -0.023    -0.022    -0.024   
                                 (0.025)   (0.025)   (0.025)  
factor(municip): 159/101         -0.051*   -0.050*   -0.051*  
                                 (0.025)   (0.025)   (0.025)  
factor(municip): 161/101         -0.090*   -0.089    -0.090*  
                                 (0.046)   (0.046)   (0.046)  
factor(municip): 163/101         -0.003    -0.002    -0.003   
                                 (0.033)   (0.033)   (0.033)  
factor(municip): 165/101         -0.179*** -0.179*** -0.181***
                                 (0.046)   (0.046)   (0.046)  
factor(municip): 167/101         -0.045    -0.045    -0.045   
                                 (0.034)   (0.034)   (0.034)  
factor(municip): 169/101         -0.052    -0.051    -0.052   
                                 (0.038)   (0.038)   (0.038)  
factor(municip): 171/101         -0.108    -0.106    -0.104   
                                 (0.081)   (0.082)   (0.082)  
factor(municip): 173/101         -0.016    -0.014    -0.014   
                                 (0.027)   (0.027)   (0.027)  
factor(municip): 175/101          0.010     0.011     0.010   
                                 (0.034)   (0.034)   (0.034)  
factor(municip): 181/101         -0.030    -0.028    -0.028   
                                 (0.042)   (0.042)   (0.042)  
factor(municip): 183/101          0.037     0.037     0.035   
                                 (0.046)   (0.046)   (0.046)  
factor(municip): 185/101         -0.057    -0.056    -0.057   
                                 (0.030)   (0.030)   (0.030)  
factor(municip): 187/101         -0.099    -0.099    -0.098   
                                 (0.052)   (0.052)   (0.052)  
factor(municip): 189/101         -0.015    -0.014    -0.010   
                                 (0.054)   (0.054)   (0.054)  
factor(municip): 190/101         -0.041    -0.040    -0.040   
                                 (0.051)   (0.051)   (0.051)  
factor(municip): 201/101          0.019     0.021     0.019   
                                 (0.048)   (0.049)   (0.049)  
factor(municip): 205/101          0.024     0.026     0.025   
                                 (0.049)   (0.050)   (0.050)  
factor(municip): 207/101         -0.091    -0.089    -0.089   
                                 (0.063)   (0.064)   (0.064)  
factor(municip): 208/101          0.017     0.018     0.018   
                                 (0.052)   (0.052)   (0.052)  
factor(municip): 209/101          0.045     0.046     0.046   
                                 (0.063)   (0.064)   (0.064)  
factor(municip): 210/101          0.067     0.068     0.068   
                                 (0.056)   (0.056)   (0.056)  
factor(municip): 211/101         -0.112*   -0.112*   -0.110*  
                                 (0.048)   (0.048)   (0.048)  
factor(municip): 213/101          0.000     0.001     0.002   
                                 (0.064)   (0.064)   (0.064)  
factor(municip): 215/101          0.038     0.039     0.038   
                                 (0.067)   (0.067)   (0.067)  
factor(municip): 217/101         -0.067*   -0.067*   -0.068*  
                                 (0.027)   (0.027)   (0.027)  
factor(municip): 219/101         -0.052    -0.050    -0.050   
                                 (0.031)   (0.032)   (0.032)  
factor(municip): 221/101          0.042     0.043     0.048   
                                 (0.076)   (0.076)   (0.076)  
factor(municip): 223/101         -0.064    -0.063    -0.062   
                                 (0.035)   (0.036)   (0.036)  
factor(municip): 225/101         -0.002    -0.001    -0.003   
                                 (0.099)   (0.099)   (0.099)  
factor(municip): 227/101          0.003     0.004     0.002   
                                 (0.058)   (0.058)   (0.058)  
factor(municip): 229/101         -0.062    -0.060    -0.062   
                                 (0.099)   (0.099)   (0.099)  
factor(municip): 230/101         -0.031    -0.030    -0.031   
                                 (0.048)   (0.048)   (0.048)  
factor(municip): 231/101         -0.015    -0.012    -0.014   
                                 (0.081)   (0.082)   (0.082)  
factor(municip): 233/101          0.014     0.015     0.016   
                                 (0.058)   (0.059)   (0.059)  
factor(municip): 235/101          0.036     0.038     0.036   
                                 (0.061)   (0.061)   (0.061)  
factor(municip): 237/101         -0.058    -0.057    -0.057   
                                 (0.056)   (0.056)   (0.056)  
factor(municip): 240/101         -0.106    -0.104    -0.105   
                                 (0.056)   (0.056)   (0.056)  
factor(municip): 250/101         -0.071    -0.070    -0.072   
                                 (0.043)   (0.043)   (0.043)  
factor(municip): 251/101          0.082     0.084     0.082   
                                 (0.099)   (0.100)   (0.099)  
factor(municip): 253/101         -0.026    -0.025    -0.023   
                                 (0.030)   (0.031)   (0.031)  
factor(municip): 255/101          0.055     0.057     0.059   
                                 (0.049)   (0.050)   (0.050)  
factor(municip): 257/101          0.047     0.049     0.047   
                                 (0.081)   (0.082)   (0.082)  
factor(municip): 259/101         -0.057    -0.056    -0.055   
                                 (0.029)   (0.030)   (0.029)  
factor(municip): 260/101         -0.060    -0.060    -0.058   
                                 (0.052)   (0.052)   (0.052)  
factor(municip): 261/101         -0.057    -0.054    -0.056   
                                 (0.140)   (0.140)   (0.140)  
factor(municip): 263/101          0.070     0.072     0.077   
                                 (0.081)   (0.082)   (0.082)  
factor(municip): 265/101         -0.028    -0.027    -0.026   
                                 (0.026)   (0.026)   (0.026)  
factor(municip): 267/101          0.064     0.066     0.063   
                                 (0.099)   (0.099)   (0.099)  
factor(municip): 269/101         -0.041    -0.039    -0.039   
                                 (0.046)   (0.046)   (0.046)  
factor(municip): 270/101         -0.025    -0.024    -0.025   
                                 (0.045)   (0.045)   (0.045)  
factor(municip): 271/101         -0.028    -0.026    -0.028   
                                 (0.081)   (0.082)   (0.082)  
factor(municip): 301/101          0.055     0.054     0.057   
                                 (0.099)   (0.099)   (0.099)  
factor(municip): 303/101          0.018     0.018     0.022   
                                 (0.115)   (0.115)   (0.114)  
factor(municip): 305/101         -0.101    -0.100    -0.103   
                                 (0.114)   (0.114)   (0.114)  
factor(municip): 306/101         -0.037    -0.037    -0.037   
                                 (0.061)   (0.061)   (0.061)  
factor(municip): 307/101         -0.157    -0.157    -0.153   
                                 (0.099)   (0.099)   (0.099)  
factor(municip): 309/101         -0.015    -0.015    -0.006   
                                 (0.099)   (0.099)   (0.099)  
factor(municip): 313/101          0.015     0.016     0.019   
                                 (0.081)   (0.081)   (0.081)  
factor(municip): 315/101          0.014     0.014     0.012   
                                 (0.046)   (0.046)   (0.046)  
factor(municip): 316/101         -0.000    -0.000    -0.002   
                                 (0.042)   (0.042)   (0.042)  
factor(municip): 317/101         -0.072    -0.071    -0.069   
                                 (0.081)   (0.081)   (0.081)  
factor(municip): 319/101         -0.101    -0.101    -0.102   
                                 (0.114)   (0.114)   (0.114)  
factor(municip): 320/101         -0.112    -0.111    -0.113   
                                 (0.064)   (0.064)   (0.064)  
factor(municip): 321/101         -0.041    -0.040    -0.042   
                                 (0.089)   (0.089)   (0.089)  
factor(municip): 323/101          0.021     0.021     0.021   
                                 (0.049)   (0.049)   (0.049)  
factor(municip): 325/101          0.058     0.058     0.058   
                                 (0.051)   (0.051)   (0.051)  
factor(municip): 326/101         -0.073    -0.073    -0.073   
                                 (0.049)   (0.050)   (0.049)  
factor(municip): 327/101         -0.156*   -0.157*   -0.155*  
                                 (0.076)   (0.076)   (0.076)  
factor(municip): 329/101         -0.058    -0.057    -0.059   
                                 (0.032)   (0.032)   (0.032)  
factor(municip): 330/101          0.009     0.009     0.009   
                                 (0.040)   (0.040)   (0.040)  
factor(municip): 331/101         -0.027    -0.027    -0.029   
                                 (0.064)   (0.064)   (0.064)  
factor(municip): 333/101         -0.043    -0.043    -0.044   
                                 (0.039)   (0.039)   (0.039)  
factor(municip): 335/101         -0.046    -0.045    -0.047   
                                 (0.064)   (0.064)   (0.064)  
factor(municip): 336/101          0.006     0.007     0.010   
                                 (0.056)   (0.056)   (0.056)  
factor(municip): 337/101          0.060     0.061     0.058   
                                 (0.197)   (0.197)   (0.197)  
factor(municip): 339/101         -0.141    -0.140    -0.135   
                                 (0.099)   (0.099)   (0.099)  
factor(municip): 340/101         -0.130*   -0.129*   -0.128*  
                                 (0.053)   (0.053)   (0.053)  
factor(municip): 341/101         -0.035    -0.035    -0.037   
                                 (0.099)   (0.099)   (0.099)  
factor(municip): 343/101         -0.072    -0.072    -0.074   
                                 (0.140)   (0.140)   (0.140)  
factor(municip): 345/101         -0.215*   -0.214*   -0.212*  
                                 (0.089)   (0.089)   (0.089)  
factor(municip): 350/101         -0.054    -0.053    -0.049   
                                 (0.054)   (0.054)   (0.054)  
factor(municip): 351/101          0.018     0.018     0.017   
                                 (0.061)   (0.061)   (0.061)  
factor(municip): 353/101          0.068     0.069     0.067   
                                 (0.197)   (0.197)   (0.197)  
factor(municip): 355/101         -0.320    -0.321    -0.300   
                                 (0.197)   (0.197)   (0.197)  
factor(municip): 357/101         -0.084    -0.083    -0.081   
                                 (0.076)   (0.076)   (0.076)  
factor(municip): 360/101         -0.033    -0.034    -0.034   
                                 (0.051)   (0.052)   (0.052)  
factor(municip): 361/101          0.052     0.054     0.050   
                                 (0.197)   (0.197)   (0.197)  
factor(municip): 363/101          0.122     0.123     0.119   
                                 (0.099)   (0.099)   (0.099)  
factor(municip): 365/101         -0.170    -0.169    -0.173   
                                 (0.114)   (0.115)   (0.114)  
factor(municip): 367/101         -0.172**  -0.173**  -0.172** 
                                 (0.058)   (0.059)   (0.059)  
factor(municip): 369/101         -0.204*** -0.204*** -0.205***
                                 (0.058)   (0.058)   (0.058)  
factor(municip): 370/101         -0.015    -0.015    -0.015   
                                 (0.034)   (0.034)   (0.034)  
factor(municip): 371/101          0.025     0.025     0.021   
                                 (0.197)   (0.197)   (0.197)  
factor(municip): 373/101         -0.014    -0.014    -0.014   
                                 (0.036)   (0.036)   (0.036)  
factor(municip): 375/101          0.055     0.055     0.053   
                                 (0.089)   (0.089)   (0.089)  
factor(municip): 376/101         -0.020    -0.020    -0.022   
                                 (0.040)   (0.040)   (0.040)  
factor(municip): 377/101         -0.013    -0.012    -0.014   
                                 (0.099)   (0.099)   (0.099)  
factor(municip): 379/101          0.001     0.000     0.006   
                                 (0.114)   (0.115)   (0.115)  
factor(municip): 383/101         -0.091    -0.093    -0.098   
                                 (0.140)   (0.140)   (0.140)  
factor(municip): 385/101          0.043     0.044     0.043   
                                 (0.140)   (0.140)   (0.140)  
factor(municip): 387/101         -0.001    -0.002    -0.004   
                                 (0.089)   (0.089)   (0.089)  
factor(municip): 389/101         -0.200*   -0.199*   -0.201*  
                                 (0.081)   (0.082)   (0.082)  
factor(municip): 390/101          0.039     0.039     0.038   
                                 (0.061)   (0.061)   (0.061)  
factor(municip): 391/101         -0.399*   -0.398*   -0.400*  
                                 (0.197)   (0.197)   (0.197)  
factor(municip): 393/101          0.047     0.049     0.048   
                                 (0.197)   (0.197)   (0.197)  
factor(municip): 395/101         -0.053    -0.052    -0.049   
                                 (0.099)   (0.099)   (0.099)  
factor(municip): 397/101         -0.083    -0.083    -0.083   
                                 (0.045)   (0.045)   (0.045)  
factor(municip): 400/101         -0.040    -0.040    -0.039   
                                 (0.041)   (0.041)   (0.041)  
factor(municip): 405/101         -0.028    -0.028    -0.030   
                                 (0.197)   (0.197)   (0.197)  
factor(municip): 407/101          0.000     0.001    -0.003   
                                 (0.081)   (0.081)   (0.081)  
factor(municip): 410/101         -0.027    -0.026    -0.025   
                                 (0.056)   (0.056)   (0.056)  
factor(municip): 420/101          0.034     0.035     0.034   
                                 (0.050)   (0.050)   (0.050)  
factor(municip): 421/101          0.042     0.041     0.040   
                                 (0.081)   (0.081)   (0.081)  
factor(municip): 423/101         -0.003    -0.003    -0.005   
                                 (0.099)   (0.099)   (0.099)  
factor(municip): 425/101         -0.035    -0.034    -0.037   
                                 (0.197)   (0.197)   (0.197)  
factor(municip): 427/101         -0.257**  -0.256**  -0.258** 
                                 (0.081)   (0.082)   (0.081)  
factor(municip): 429/101         -0.288**  -0.287**  -0.284** 
                                 (0.099)   (0.099)   (0.099)  
factor(municip): 430/101         -0.110*   -0.109*   -0.110*  
                                 (0.047)   (0.047)   (0.047)  
factor(municip): 431/101         -0.112    -0.112    -0.109   
                                 (0.061)   (0.061)   (0.061)  
factor(municip): 433/101          0.061     0.062     0.060   
                                 (0.099)   (0.099)   (0.099)  
factor(municip): 435/101          0.111     0.112     0.117   
                                 (0.115)   (0.115)   (0.115)  
factor(municip): 437/101          0.058     0.059     0.061   
                                 (0.082)   (0.082)   (0.082)  
factor(municip): 439/101         -0.160    -0.158    -0.160   
                                 (0.114)   (0.115)   (0.115)  
factor(municip): 440/101         -0.047    -0.047    -0.047   
                                 (0.058)   (0.058)   (0.058)  
factor(municip): 441/101         -0.100    -0.098    -0.095   
                                 (0.099)   (0.099)   (0.099)  
factor(municip): 443/101         -0.203    -0.202    -0.182   
                                 (0.197)   (0.197)   (0.198)  
factor(municip): 445/101          0.049     0.050     0.050   
                                 (0.052)   (0.053)   (0.053)  
factor(municip): 447/101         -0.119    -0.119    -0.121   
                                 (0.076)   (0.076)   (0.076)  
factor(municip): 449/101          0.050     0.050     0.051   
                                 (0.051)   (0.051)   (0.051)  
factor(municip): 450/101         -0.040    -0.040    -0.041   
                                 (0.045)   (0.045)   (0.045)  
factor(municip): 451/101         -0.035    -0.033    -0.030   
                                 (0.089)   (0.089)   (0.089)  
factor(municip): 461/101         -0.013    -0.013    -0.014   
                                 (0.017)   (0.017)   (0.017)  
factor(municip): 471/101         -0.041    -0.040    -0.038   
                                 (0.089)   (0.089)   (0.089)  
factor(municip): 473/101         -0.026    -0.026    -0.028   
                                 (0.064)   (0.064)   (0.064)  
factor(municip): 475/101         -0.123    -0.124    -0.126   
                                 (0.081)   (0.081)   (0.081)  
factor(municip): 477/101         -0.036    -0.036    -0.034   
                                 (0.081)   (0.081)   (0.081)  
factor(municip): 479/101         -0.001    -0.001    -0.002   
                                 (0.028)   (0.028)   (0.028)  
factor(municip): 480/101         -0.031    -0.030    -0.028   
                                 (0.071)   (0.071)   (0.071)  
factor(municip): 481/101         -0.059    -0.060    -0.061   
                                 (0.197)   (0.197)   (0.197)  
factor(municip): 482/101         -0.014    -0.014    -0.016   
                                 (0.099)   (0.099)   (0.099)  
factor(municip): 483/101         -0.014    -0.013    -0.015   
                                 (0.089)   (0.089)   (0.089)  
factor(municip): 485/101         -0.014    -0.013    -0.014   
                                 (0.099)   (0.099)   (0.099)  
factor(municip): 487/101         -0.109    -0.109    -0.112   
                                 (0.197)   (0.197)   (0.197)  
factor(municip): 489/101         -0.011    -0.009    -0.012   
                                 (0.089)   (0.089)   (0.089)  
factor(municip): 491/101          0.032     0.032     0.031   
                                 (0.099)   (0.099)   (0.099)  
factor(municip): 492/101         -0.093    -0.092    -0.091   
                                 (0.071)   (0.071)   (0.071)  
factor(municip): 495/101          0.016     0.018     0.015   
                                 (0.114)   (0.114)   (0.114)  
factor(municip): 497/101         -0.024    -0.022    -0.025   
                                 (0.081)   (0.082)   (0.082)  
factor(municip): 499/101          0.071     0.072     0.070   
                                 (0.114)   (0.114)   (0.114)  
factor(municip): 501/101         -0.009    -0.007    -0.010   
                                 (0.115)   (0.115)   (0.115)  
factor(municip): 503/101         -0.061    -0.061    -0.057   
                                 (0.076)   (0.076)   (0.076)  
factor(municip): 505/101         -0.023    -0.022    -0.025   
                                 (0.197)   (0.197)   (0.197)  
factor(municip): 507/101         -0.124    -0.123    -0.124   
                                 (0.099)   (0.100)   (0.100)  
factor(municip): 509/101         -0.009    -0.008    -0.009   
                                 (0.099)   (0.099)   (0.099)  
factor(municip): 510/101          0.019     0.020     0.018   
                                 (0.039)   (0.039)   (0.039)  
factor(municip): 511/101         -0.135    -0.134    -0.131   
                                 (0.099)   (0.099)   (0.099)  
factor(municip): 513/101         -0.019    -0.018    -0.022   
                                 (0.099)   (0.099)   (0.099)  
factor(municip): 515/101          0.026     0.026     0.026   
                                 (0.041)   (0.041)   (0.041)  
factor(municip): 517/101          0.009     0.009     0.007   
                                 (0.140)   (0.140)   (0.140)  
factor(municip): 519/101         -0.031    -0.030    -0.032   
                                 (0.082)   (0.082)   (0.082)  
factor(municip): 521/101          0.052     0.053     0.052   
                                 (0.114)   (0.115)   (0.114)  
factor(municip): 523/101         -0.051    -0.050    -0.052   
                                 (0.071)   (0.071)   (0.071)  
factor(municip): 525/101         -0.158    -0.158    -0.149   
                                 (0.140)   (0.140)   (0.140)  
factor(municip): 527/101         -0.116    -0.114    -0.116   
                                 (0.114)   (0.115)   (0.115)  
factor(municip): 529/101         -0.020    -0.019    -0.019   
                                 (0.064)   (0.064)   (0.064)  
factor(municip): 530/101          0.012     0.013     0.014   
                                 (0.054)   (0.054)   (0.054)  
factor(municip): 531/101         -0.003    -0.002    -0.004   
                                 (0.114)   (0.114)   (0.114)  
factor(municip): 533/101          0.022     0.024     0.022   
                                 (0.114)   (0.115)   (0.115)  
factor(municip): 535/101          0.010     0.012     0.011   
                                 (0.114)   (0.115)   (0.115)  
factor(municip): 537/101         -0.020    -0.020    -0.022   
                                 (0.045)   (0.045)   (0.045)  
factor(municip): 539/101         -0.032    -0.031    -0.033   
                                 (0.082)   (0.082)   (0.082)  
factor(municip): 540/101          0.032     0.033     0.033   
                                 (0.035)   (0.035)   (0.035)  
factor(municip): 541/101          0.010     0.009     0.008   
                                 (0.071)   (0.071)   (0.071)  
factor(municip): 543/101         -0.044    -0.043    -0.045   
                                 (0.054)   (0.054)   (0.054)  
factor(municip): 545/101         -0.027    -0.026    -0.027   
                                 (0.051)   (0.051)   (0.051)  
factor(municip): 550/101          0.026     0.027     0.025   
                                 (0.053)   (0.053)   (0.053)  
factor(municip): 551/101         -0.206*   -0.205*   -0.201*  
                                 (0.099)   (0.099)   (0.099)  
factor(municip): 553/101         -0.072    -0.070    -0.072   
                                 (0.140)   (0.140)   (0.140)  
factor(municip): 555/101          0.150     0.151     0.157   
                                 (0.114)   (0.114)   (0.114)  
factor(municip): 557/101         -0.103    -0.102    -0.101   
                                 (0.076)   (0.076)   (0.076)  
factor(municip): 559/101          0.163     0.164     0.162   
                                 (0.114)   (0.114)   (0.114)  
factor(municip): 561/101         -0.009    -0.009    -0.009   
                                 (0.021)   (0.021)   (0.021)  
factor(municip): 563/101         -0.049    -0.048    -0.045   
                                 (0.056)   (0.056)   (0.056)  
factor(municip): 565/101         -0.072    -0.071    -0.070   
                                 (0.071)   (0.071)   (0.071)  
factor(municip): 567/101          0.058     0.060     0.058   
                                 (0.115)   (0.115)   (0.115)  
factor(municip): 569/101          0.065     0.066     0.069   
                                 (0.089)   (0.089)   (0.089)  
factor(municip): 571/101         -0.127    -0.126    -0.125   
                                 (0.067)   (0.067)   (0.067)  
factor(municip): 573/101         -0.065    -0.064    -0.064   
                                 (0.040)   (0.040)   (0.040)  
factor(municip): 575/101         -0.017    -0.016    -0.015   
                                 (0.033)   (0.033)   (0.033)  
factor(municip): 577/101         -0.033    -0.031    -0.033   
                                 (0.076)   (0.076)   (0.076)  
factor(municip): 580/101         -0.015    -0.014    -0.012   
                                 (0.040)   (0.040)   (0.040)  
factor(municip): 601/101         -0.250**  -0.249**  -0.246** 
                                 (0.089)   (0.089)   (0.089)  
factor(municip): 603/101          0.051     0.053     0.057   
                                 (0.099)   (0.100)   (0.100)  
factor(municip): 605/101          0.003     0.005     0.005   
                                 (0.056)   (0.057)   (0.057)  
factor(municip): 607/101          0.015     0.016     0.015   
                                 (0.026)   (0.026)   (0.026)  
factor(municip): 609/101         -0.052    -0.050    -0.048   
                                 (0.089)   (0.089)   (0.089)  
factor(municip): 611/101          0.042     0.043     0.046   
                                 (0.089)   (0.089)   (0.089)  
factor(municip): 613/101         -0.041    -0.039    -0.041   
                                 (0.058)   (0.059)   (0.059)  
factor(municip): 615/101         -0.018    -0.018    -0.018   
                                 (0.023)   (0.023)   (0.023)  
factor(municip): 617/101          0.085     0.087     0.093   
                                 (0.114)   (0.115)   (0.115)  
factor(municip): 619/101         -0.074    -0.073    -0.075   
                                 (0.071)   (0.071)   (0.071)  
factor(municip): 621/101         -0.048*   -0.047*   -0.047*  
                                 (0.023)   (0.023)   (0.023)  
factor(municip): 623/101         -0.167    -0.165    -0.167   
                                 (0.089)   (0.089)   (0.089)  
factor(municip): 625/101         -0.053    -0.051    -0.053   
                                 (0.197)   (0.197)   (0.197)  
factor(municip): 627/101         -0.030    -0.029    -0.026   
                                 (0.099)   (0.099)   (0.099)  
factor(municip): 629/101          0.003     0.005     0.002   
                                 (0.089)   (0.089)   (0.089)  
factor(municip): 630/101         -0.058    -0.057    -0.058   
                                 (0.036)   (0.036)   (0.036)  
factor(municip): 631/101         -0.035    -0.034    -0.035   
                                 (0.035)   (0.035)   (0.035)  
factor(municip): 651/101          0.134     0.135     0.138   
                                 (0.099)   (0.100)   (0.100)  
factor(municip): 653/101         -0.041    -0.039    -0.035   
                                 (0.076)   (0.076)   (0.076)  
factor(municip): 655/101          0.015     0.016     0.026   
                                 (0.140)   (0.140)   (0.140)  
factor(municip): 657/101         -0.009    -0.009    -0.009   
                                 (0.024)   (0.024)   (0.024)  
factor(municip): 659/101         -0.047    -0.045    -0.047   
                                 (0.140)   (0.140)   (0.140)  
factor(municip): 661/101         -0.059*   -0.058*   -0.059*  
                                 (0.029)   (0.029)   (0.029)  
factor(municip): 663/101         -0.016    -0.015    -0.014   
                                 (0.053)   (0.053)   (0.053)  
factor(municip): 665/101         -0.024    -0.023    -0.022   
                                 (0.041)   (0.041)   (0.041)  
factor(municip): 667/101          0.071     0.073     0.071   
                                 (0.076)   (0.076)   (0.076)  
factor(municip): 669/101         -0.014    -0.012    -0.014   
                                 (0.099)   (0.099)   (0.099)  
factor(municip): 671/101         -0.014    -0.014    -0.016   
                                 (0.044)   (0.044)   (0.044)  
factor(municip): 673/101         -0.236*   -0.235*   -0.238*  
                                 (0.114)   (0.115)   (0.115)  
factor(municip): 675/101          0.132     0.134     0.144   
                                 (0.140)   (0.140)   (0.140)  
factor(municip): 677/101          0.079     0.080     0.078   
                                 (0.099)   (0.099)   (0.099)  
factor(municip): 679/101          0.035     0.037     0.048   
                                 (0.089)   (0.089)   (0.089)  
factor(municip): 681/101          0.012     0.014     0.011   
                                 (0.082)   (0.082)   (0.082)  
factor(municip): 683/101          0.096     0.097     0.103   
                                 (0.115)   (0.115)   (0.115)  
factor(municip): 685/101          0.156     0.157     0.155   
                                 (0.140)   (0.140)   (0.140)  
factor(municip): 701/101          0.025     0.025     0.024   
                                 (0.076)   (0.076)   (0.076)  
factor(municip): 703/101          0.087     0.089     0.089   
                                 (0.064)   (0.064)   (0.064)  
factor(municip): 706/101          0.029     0.030     0.031   
                                 (0.041)   (0.041)   (0.041)  
factor(municip): 707/101          0.012     0.012     0.011   
                                 (0.043)   (0.043)   (0.043)  
factor(municip): 709/101         -0.064    -0.063    -0.065   
                                 (0.082)   (0.082)   (0.082)  
factor(municip): 710/101         -0.013    -0.011    -0.012   
                                 (0.038)   (0.038)   (0.038)  
factor(municip): 711/101         -0.024    -0.022    -0.022   
                                 (0.071)   (0.071)   (0.071)  
factor(municip): 713/101         -0.138    -0.136    -0.138   
                                 (0.071)   (0.071)   (0.071)  
factor(municip): 715/101         -0.101    -0.100    -0.102   
                                 (0.064)   (0.064)   (0.064)  
factor(municip): 717/101          0.036     0.037     0.039   
                                 (0.081)   (0.082)   (0.082)  
factor(municip): 719/101         -0.051    -0.049    -0.041   
                                 (0.140)   (0.140)   (0.140)  
factor(municip): 721/101          0.004     0.005     0.003   
                                 (0.089)   (0.089)   (0.089)  
factor(municip): 723/101          0.027     0.028     0.026   
                                 (0.064)   (0.064)   (0.064)  
factor(municip): 725/101         -0.073    -0.073    -0.064   
                                 (0.140)   (0.140)   (0.140)  
factor(municip): 727/101         -0.004    -0.003    -0.005   
                                 (0.044)   (0.044)   (0.044)  
factor(municip): 729/101          0.008     0.009     0.007   
                                 (0.081)   (0.082)   (0.082)  
factor(municip): 730/101         -0.046    -0.045    -0.046   
                                 (0.029)   (0.029)   (0.029)  
factor(municip): 731/101         -0.025    -0.025    -0.023   
                                 (0.029)   (0.029)   (0.029)  
factor(municip): 733/101         -0.003    -0.002    -0.003   
                                 (0.081)   (0.082)   (0.082)  
factor(municip): 735/101         -0.034    -0.034    -0.036   
                                 (0.114)   (0.115)   (0.114)  
factor(municip): 737/101          0.019     0.021     0.021   
                                 (0.058)   (0.059)   (0.059)  
factor(municip): 739/101          0.003     0.005     0.003   
                                 (0.071)   (0.071)   (0.071)  
factor(municip): 740/101          0.033     0.034     0.033   
                                 (0.029)   (0.030)   (0.030)  
factor(municip): 741/101         -0.123    -0.122    -0.124   
                                 (0.197)   (0.197)   (0.197)  
factor(municip): 743/101         -0.038    -0.037    -0.036   
                                 (0.030)   (0.030)   (0.030)  
factor(municip): 745/101         -0.105    -0.104    -0.106   
                                 (0.058)   (0.058)   (0.058)  
factor(municip): 746/101         -0.001     0.001     0.001   
                                 (0.038)   (0.039)   (0.039)  
factor(municip): 747/101          0.027     0.028     0.027   
                                 (0.081)   (0.082)   (0.082)  
factor(municip): 749/101         -0.089    -0.088    -0.084   
                                 (0.099)   (0.099)   (0.099)  
factor(municip): 751/101         -0.014    -0.014    -0.014   
                                 (0.015)   (0.015)   (0.015)  
factor(municip): 756/101         -0.087    -0.086    -0.085   
                                 (0.051)   (0.051)   (0.051)  
factor(municip): 760/101          0.001     0.002     0.003   
                                 (0.040)   (0.040)   (0.040)  
factor(municip): 761/101          0.015     0.017     0.025   
                                 (0.140)   (0.140)   (0.140)  
factor(municip): 763/101         -0.075    -0.074    -0.076   
                                 (0.089)   (0.089)   (0.089)  
factor(municip): 765/101         -0.146    -0.145    -0.135   
                                 (0.140)   (0.140)   (0.140)  
factor(municip): 766/101          0.043     0.045     0.045   
                                 (0.046)   (0.046)   (0.046)  
factor(municip): 767/101          0.131     0.133     0.136   
                                 (0.089)   (0.089)   (0.089)  
factor(municip): 769/101         -0.090    -0.089    -0.091   
                                 (0.197)   (0.197)   (0.197)  
factor(municip): 771/101          0.017     0.018     0.017   
                                 (0.076)   (0.076)   (0.076)  
factor(municip): 773/101         -0.082    -0.082    -0.081   
                                 (0.045)   (0.045)   (0.045)  
factor(municip): 775/101          0.163     0.165     0.165   
                                 (0.197)   (0.197)   (0.197)  
factor(municip): 777/101         -0.178*   -0.177*   -0.179*  
                                 (0.089)   (0.089)   (0.089)  
factor(municip): 779/101         -0.054    -0.053    -0.054   
                                 (0.032)   (0.032)   (0.032)  
factor(municip): 781/101          0.071     0.072     0.069   
                                 (0.140)   (0.140)   (0.140)  
factor(municip): 783/101         -0.015    -0.013    -0.014   
                                 (0.114)   (0.115)   (0.115)  
factor(municip): 785/101          0.024     0.025     0.026   
                                 (0.071)   (0.071)   (0.071)  
factor(municip): 787/101         -0.034    -0.034    -0.033   
                                 (0.033)   (0.033)   (0.033)  
factor(municip): 789/101          0.075     0.075     0.094   
                                 (0.197)   (0.197)   (0.198)  
factor(municip): 791/101          0.003     0.004     0.004   
                                 (0.027)   (0.027)   (0.027)  
factor(municip): 793/101         -0.046    -0.045    -0.047   
                                 (0.089)   (0.089)   (0.089)  
factor(municip): 801/101         -0.052    -0.050    -0.053   
                                 (0.140)   (0.140)   (0.140)  
factor(municip): 803/101          0.005     0.006     0.004   
                                 (0.197)   (0.197)   (0.197)  
factor(municip): 805/101         -0.018    -0.017    -0.016   
                                 (0.047)   (0.047)   (0.047)  
factor(municip): 807/101         -0.063    -0.062    -0.062   
                                 (0.064)   (0.064)   (0.064)  
factor(municip): 809/101          0.192     0.194     0.192   
                                 (0.197)   (0.197)   (0.197)  
factor(municip): 810/101         -0.148*** -0.147*** -0.148***
                                 (0.044)   (0.044)   (0.044)  
factor(municip): 811/101         -0.109    -0.108    -0.110   
                                 (0.064)   (0.064)   (0.064)  
factor(municip): 813/101         -0.015    -0.014    -0.014   
                                 (0.028)   (0.028)   (0.028)  
factor(municip): 815/101         -0.069    -0.068    -0.071   
                                 (0.099)   (0.099)   (0.099)  
factor(municip): 817/101         -0.024    -0.023    -0.025   
                                 (0.114)   (0.114)   (0.114)  
factor(municip): 819/101          0.040     0.041     0.039   
                                 (0.056)   (0.056)   (0.056)  
factor(municip): 820/101          0.001     0.002     0.002   
                                 (0.061)   (0.061)   (0.061)  
factor(municip): 821/101          0.008     0.008     0.007   
                                 (0.041)   (0.041)   (0.041)  
factor(municip): 823/101         -0.019    -0.018    -0.020   
                                 (0.099)   (0.099)   (0.099)  
factor(municip): 827/101         -0.314**  -0.314**  -0.316** 
                                 (0.099)   (0.099)   (0.099)  
factor(municip): 829/101         -0.119    -0.119    -0.121   
                                 (0.099)   (0.099)   (0.099)  
factor(municip): 831/101         -0.013    -0.012    -0.013   
                                 (0.081)   (0.082)   (0.082)  
factor(municip): 833/101         -0.116    -0.115    -0.117   
                                 (0.115)   (0.115)   (0.115)  
factor(municip): 835/101         -0.062    -0.062    -0.063   
                                 (0.114)   (0.114)   (0.114)  
factor(municip): 837/101          0.023     0.024     0.023   
                                 (0.071)   (0.071)   (0.071)  
factor(municip): 839/101          0.082     0.084     0.084   
                                 (0.082)   (0.082)   (0.082)  
factor(municip): 840/101         -0.042    -0.040    -0.039   
                                 (0.054)   (0.055)   (0.055)  
factor(municip): 841/101         -0.005    -0.004     0.000   
                                 (0.061)   (0.061)   (0.061)  
factor(municip): 843/101         -0.057    -0.055    -0.045   
                                 (0.140)   (0.140)   (0.140)  
factor(municip): 845/101         -0.085    -0.083    -0.079   
                                 (0.099)   (0.100)   (0.100)  
factor(municip): 846/101         -0.026    -0.026    -0.024   
                                 (0.050)   (0.050)   (0.050)  
factor(municip): 847/101          0.042     0.043     0.041   
                                 (0.056)   (0.056)   (0.056)  
factor(municip): 849/101         -0.045    -0.044    -0.045   
                                 (0.036)   (0.036)   (0.036)  
factor(municip): 851/101         -0.016    -0.015    -0.015   
                                 (0.017)   (0.017)   (0.017)  
factor(municip): 860/101         -0.002    -0.001    -0.002   
                                 (0.033)   (0.033)   (0.033)  
factor(municip): 861/101         -0.020    -0.018    -0.018   
                                 (0.061)   (0.061)   (0.061)  
factor(essround): 2/1             0.015     0.014     0.014   
                                 (0.009)   (0.009)   (0.009)  
factor(essround): 3/1             0.067***  0.067***  0.067***
                                 (0.008)   (0.009)   (0.009)  
factor(essround): 4/1            -0.077*** -0.077*** -0.077***
                                 (0.009)   (0.010)   (0.010)  
factor(essround): 5/1            -0.150*** -0.150*** -0.150***
                                 (0.009)   (0.009)   (0.009)  
Population (80m)                 -0.000**  -0.000**  -0.000** 
                                 (0.000)   (0.000)   (0.000)  
Unemployed                       -0.052*** -0.052*** -0.050***
                                 (0.009)   (0.009)   (0.009)  
Age                               0.000     0.000     0.000   
                                 (0.000)   (0.000)   (0.000)  
Years of education                0.002*    0.002*    0.002*  
                                 (0.001)   (0.001)   (0.001)  
Male                              0.045***  0.045***  0.046***
                                 (0.005)   (0.005)   (0.005)  
Income (in 100k DKR)              0.004*    0.004*    0.004*  
                                 (0.002)   (0.002)   (0.002)  
Years lived at location           0.000     0.000     0.000   
                                 (0.000)   (0.000)   (0.000)  
Right--wing                       0.166***  0.166***  0.166***
                                 (0.013)   (0.013)   (0.013)  
Newspaper exposure               -0.008    -0.008    -0.009   
                                 (0.017)   (0.017)   (0.017)  
Political interest               -0.009*   -0.009*   -0.009*  
                                 (0.004)   (0.004)   (0.004)  
Local Unemployment Rate (80m)    -0.070*   -0.072*   -0.064*  
                                 (0.030)   (0.031)   (0.031)  
Local Unemployment Rate (2500m)             0.029     0.023   
                                           (0.103)   (0.103)  
Household Unemployment                               -0.023*  
                                                     (0.010)  
--------------------------------------------------------------
R-squared                            0.223     0.223     0.223
adj. R-squared                       0.179     0.178     0.179
sigma                                0.197     0.197     0.197
F                                    5.056     5.040     5.042
p                                    0.000     0.000     0.000
Log-likelihood                    1435.195  1435.238  1437.776
Deviance                           223.152   223.149   222.963
AIC                              -2212.389 -2210.476 -2213.551
BIC                                 -2.757     5.872     9.514
N                                 6101      6101      6101    
==============================================================

> ######################################
> # Results and code for Figure 2	#
> # Sensitivity to size of context	#
> 
> 
> ueb.names<-paste("ueb",c(80,130,180,250,350,500,750,1000,1250,1500,1750,2000,2250,2500),sep="")

> popdens.names<-paste("popdens",c(80,130,180,250,350,500,750,1000,1250,1500,1750,2000,2250,2500),sep="")

> eff.out2<-matrix(NA,nr=14,nc=2)

> eff.out1<-matrix(NA,nr=14,nc=2)

> for (i in 1:14){
+ 
+ lmfit2<-lm(substitute(stfeco~sub1+sub2+factor(municip)+factor(essround)+unemplo4+ageR+eduresu+maleR+pincome2000_1k+dayslived_y+lrscale+nwsptot+polintr+partunemplo4
+ 	,list(sub1=as.name(ueb.names[i]),sub2=as.name(popdens.names[i])))
+ 	,data=dat[dat$popdens80>9,])
+ 
+ eff.out2[i,1]<-coef(lmfit2)[ueb.names[i]]
+ eff.out2[i,2]<-sqrt(diag(vcov(lmfit2)))[ueb.names[i]]
+ 
+ lmfit1<-lm(substitute(stfeco~sub1+sub2+factor(essround)+unemplo4+ageR+eduresu+maleR+pincome2000_1k+dayslived_y+lrscale+nwsptot+polintr+partunemplo4
+ 	,list(sub1=as.name(ueb.names[i]),sub2=as.name(popdens.names[i])))
+ 	,data=dat[dat$popdens80>9,])
+ 
+ eff.out1[i,1]<-coef(lmfit1)[ueb.names[i]]
+ eff.out1[i,2]<-sqrt(diag(vcov(lmfit1)))[ueb.names[i]]
+ 
+ }

> eff.out2
              [,1]       [,2]
 [1,] -0.062748600 0.03043477
 [2,] -0.107061101 0.03571154
 [3,] -0.099242120 0.03883389
 [4,] -0.094657469 0.04222879
 [5,] -0.085935423 0.04629927
 [6,] -0.069279218 0.05120137
 [7,] -0.072095340 0.05868277
 [8,] -0.064264109 0.06539287
 [9,] -0.035218735 0.07142424
[10,] -0.006186597 0.07676570
[11,] -0.006037732 0.08203864
[12,]  0.006669111 0.08852484
[13,] -0.021474567 0.09589039
[14,] -0.016511508 0.10247187

> eff.out1
             [,1]       [,2]
 [1,] -0.07897674 0.02890104
 [2,] -0.11823016 0.03376806
 [3,] -0.11805768 0.03666906
 [4,] -0.11336919 0.03976290
 [5,] -0.10483421 0.04333281
 [6,] -0.09129040 0.04743394
 [7,] -0.09308730 0.05329363
 [8,] -0.09037972 0.05843925
 [9,] -0.07298385 0.06269749
[10,] -0.06045364 0.06616864
[11,] -0.06912821 0.06945004
[12,] -0.07363016 0.07336704
[13,] -0.10290908 0.07764632
[14,] -0.10833752 0.08129965

> #weighting estimates with 1st to 99th percentile variation.
> w<-matrix(NA,nr=length(ueb.names),nc=2)

> for (i in 1:length(ueb.names))w[i,]<-quantile(dat[dat$popdens80>9,ueb.names[i]],c(.01,.99),na.rm=T)

> w.out<-apply(w,1,diff)

> (eff.out2<-eff.out2*w.out)
              [,1]       [,2]
 [1,] -0.027339974 0.01326063
 [2,] -0.044278055 0.01476949
 [3,] -0.037938722 0.01484559
 [4,] -0.031862741 0.01421467
 [5,] -0.027181904 0.01464475
 [6,] -0.019906809 0.01471229
 [7,] -0.017925985 0.01459105
 [8,] -0.014178761 0.01442780
 [9,] -0.007235623 0.01467398
[10,] -0.001225502 0.01520651
[11,] -0.001168284 0.01587424
[12,]  0.001247571 0.01656009
[13,] -0.003763415 0.01680478
[14,] -0.002745808 0.01704072

> (eff.out1<-eff.out1*w.out)
             [,1]       [,2]
 [1,] -0.03441068 0.01259237
 [2,] -0.04889733 0.01396571
 [3,] -0.04513162 0.01401801
 [4,] -0.03816131 0.01338463
 [5,] -0.03315971 0.01370643
 [6,] -0.02623154 0.01362975
 [7,] -0.02314548 0.01325108
 [8,] -0.01994072 0.01289361
 [9,] -0.01499440 0.01288108
[10,] -0.01197525 0.01310734
[11,] -0.01337611 0.01343838
[12,] -0.01377378 0.01372456
[13,] -0.01803480 0.01360751
[14,] -0.01801616 0.01351986

> ######################################
> #	PLOTTING				#
> #	With/without MUNICI FEs		#
> 
> #pdf('ueb-agg-weighted.pdf',width=6.5,height=4)
> setEPS()

> postscript("ueb-agg-weighted.eps",width=6.5,height=4)

> jit<-c(rep(10,5),15,rep(20,8))

> xax<-c(80,130,180,250,350,500,750,1000,1250,1500,1750,2000,2250,2500)

> par(mfrow=c(1,1),omd=c(0.05,1,0,1))

> plot(xax,eff.out2[,1],main="",yaxt="n",col="white",
+ 	ylim=c(-.1,.05),xaxt="n",
+ 	ylab="Marginal Effect of Local Unemployment Rate",
+ 	xlab="Level of aggregation of contextual predictor (radius in meters)",
+ 	cex.lab=1)

> segments(xax-jit,eff.out2[,1],xax-jit,eff.out2[,1]+1.97*eff.out2[,2],lwd=1)

> segments(xax-jit,eff.out2[,1],xax-jit,eff.out2[,1]-1.97*eff.out2[,2],lwd=1)

> abline(h=0,lty=2)

> points(xax-jit,eff.out2[,1],lwd=1,pch=21,bg="black",col="black",cex=1.25)

> axis(2,las=2,at=round(seq(-.1,.05,.05),2),cex.axis=1)

> legend("bottomleft",lwd=1.5,lty=1,bty="n",legend=c("95% CI"),cex=1)

> legend("bottomright",pch=21,pt.bg=c("black","white"),legend=c("With municipality fixed effects","Without municipality fixed effects"),bty="n",
+ 	cex=1,pt.cex=1,col="black")

> #ADDING W/O MUNICI FEs
> segments(xax+jit,eff.out1[,1],xax+jit,eff.out1[,1]+1.97*eff.out1[,2],lwd=1,col="black")

> segments(xax+jit,eff.out1[,1],xax+jit,eff.out1[,1]-1.97*eff.out1[,2],lwd=1,col="black")

> abline(h=0,lty=2)

> points(xax+jit,eff.out1[,1],lwd=1,pch=21,bg="white",col="black",cex=1.25)

> pos<-c(0,250,500,750,1000,1250,1500,1750,2000,2250,2500)

> axis(1,at=pos,labels=NA)

> text(pos[-1],par()$usr[3]-.06*(par()$usr[4]-par()$usr[3]),labels=c("250m","500m","750m","1000m","1250m","1500m","1750m","2000m","2250m","2500m"),
+ 	adj=1,srt=45,xpd=T,cex=1)

> dev.off()
null device 
          1 

> ##################################################
> # Results and code for Table 2 and Figure 3	#
> 
> fit2FE<-lm(stfeco~factor(municip)+factor(essround)+unemplo4+ageR+eduresu+maleR+lrscale+polintr+oneunemplo4
+ 	+pincome2000_1k+dayslived_y+popdens80+ueb80*nwsptot+munemploy*nwsptot
+ 	,data=dat[dat$popdens80>9,])

> fit1FE<-lm(stfeco~factor(municip)+factor(essround)+unemplo4+ageR+eduresu+maleR+lrscale+polintr+oneunemplo4
+ 	+pincome2000_1k+dayslived_y+popdens80+ueb80*nwsptot+popdens2500+ueb2500*nwsptot
+ 	,data=dat[rownames(dat)%in%rownames(fit2FE$model),])

> tab<-mtable(fit1FE,fit2FE)

> tab<-relabel(tab,
+ 	"ueb80"="Local Unemployment Rate (80m)",
+ 	ueb2500="Local Unemployment Rate (2500m)",
+ 	ueb2500="Local Unemployment Rate (2500m)",
+ 	oneunemplo4="Household Unemployment",
+ 	lrscale="Right--wing",
+ 	polintr="Political interest",
+ 	nwsptot="Newspaper exposure",
+ 	dayslived_y="Years lived at location",
+ 	pincome2000_1k="Income (in 100k DKR)",
+ 	ageR="Age",
+ 	unemplo4="Unemployed",
+ 	popdens80="Population (80m)",
+ 	eduresu="Years of education",
+ 	"maleR: Male/Female"="Male"
+ 	)

> tab

Calls:
fit1FE: lm(formula = stfeco ~ factor(municip) + factor(essround) + unemplo4 + 
    ageR + eduresu + maleR + lrscale + polintr + oneunemplo4 + 
    pincome2000_1k + dayslived_y + popdens80 + ueb80 * nwsptot + 
    popdens2500 + ueb2500 * nwsptot, data = dat[rownames(dat) %in% 
    rownames(fit2FE$model), ])
fit2FE: lm(formula = stfeco ~ factor(municip) + factor(essround) + unemplo4 + 
    ageR + eduresu + maleR + lrscale + polintr + oneunemplo4 + 
    pincome2000_1k + dayslived_y + popdens80 + ueb80 * nwsptot + 
    munemploy * nwsptot, data = dat[dat$popdens80 > 9, ])

====================================================
                                  fit1FE    fit2FE  
----------------------------------------------------
(Intercept)                       0.594***  0.620***
                                 (0.033)   (0.038)  
factor(municip): 147/101          0.013    -0.005   
                                 (0.023)   (0.022)  
factor(municip): 151/101         -0.086*   -0.087*  
                                 (0.039)   (0.038)  
factor(municip): 153/101         -0.008    -0.002   
                                 (0.041)   (0.039)  
factor(municip): 155/101         -0.031    -0.035   
                                 (0.046)   (0.046)  
factor(municip): 157/101         -0.032    -0.042   
                                 (0.026)   (0.028)  
factor(municip): 159/101         -0.059*   -0.061*  
                                 (0.026)   (0.026)  
factor(municip): 161/101         -0.100*   -0.098*  
                                 (0.047)   (0.046)  
factor(municip): 163/101         -0.013    -0.012   
                                 (0.035)   (0.034)  
factor(municip): 165/101         -0.194*** -0.184***
                                 (0.048)   (0.046)  
factor(municip): 167/101         -0.053    -0.051   
                                 (0.035)   (0.034)  
factor(municip): 169/101         -0.063    -0.058   
                                 (0.040)   (0.038)  
factor(municip): 171/101         -0.113    -0.123   
                                 (0.083)   (0.083)  
factor(municip): 173/101         -0.026    -0.033   
                                 (0.029)   (0.030)  
factor(municip): 175/101          0.004     0.002   
                                 (0.035)   (0.035)  
factor(municip): 181/101         -0.052    -0.058   
                                 (0.048)   (0.048)  
factor(municip): 183/101          0.023     0.037   
                                 (0.048)   (0.046)  
factor(municip): 185/101         -0.068*   -0.065*  
                                 (0.032)   (0.030)  
factor(municip): 187/101         -0.108*   -0.114*  
                                 (0.054)   (0.054)  
factor(municip): 189/101         -0.023    -0.031   
                                 (0.056)   (0.056)  
factor(municip): 190/101         -0.052    -0.051   
                                 (0.052)   (0.052)  
factor(municip): 201/101          0.004    -0.003   
                                 (0.050)   (0.050)  
factor(municip): 205/101          0.026     0.018   
                                 (0.062)   (0.063)  
factor(municip): 207/101         -0.101    -0.105   
                                 (0.065)   (0.065)  
factor(municip): 208/101         -0.009    -0.009   
                                 (0.060)   (0.059)  
factor(municip): 209/101          0.015     0.015   
                                 (0.068)   (0.068)  
factor(municip): 210/101          0.051     0.049   
                                 (0.058)   (0.057)  
factor(municip): 211/101         -0.067    -0.060   
                                 (0.065)   (0.064)  
factor(municip): 213/101         -0.033    -0.030   
                                 (0.069)   (0.068)  
factor(municip): 215/101          0.022     0.024   
                                 (0.068)   (0.068)  
factor(municip): 217/101         -0.079**  -0.077** 
                                 (0.030)   (0.028)  
factor(municip): 219/101         -0.062    -0.065   
                                 (0.035)   (0.034)  
factor(municip): 221/101          0.036     0.036   
                                 (0.077)   (0.076)  
factor(municip): 223/101         -0.075*   -0.082*  
                                 (0.038)   (0.039)  
factor(municip): 225/101         -0.061    -0.059   
                                 (0.115)   (0.115)  
factor(municip): 227/101         -0.011    -0.013   
                                 (0.060)   (0.060)  
factor(municip): 229/101         -0.200    -0.198   
                                 (0.141)   (0.140)  
factor(municip): 230/101         -0.047    -0.051   
                                 (0.050)   (0.050)  
factor(municip): 231/101         -0.069    -0.079   
                                 (0.091)   (0.091)  
factor(municip): 233/101          0.002    -0.007   
                                 (0.060)   (0.061)  
factor(municip): 235/101          0.021     0.013   
                                 (0.083)   (0.083)  
factor(municip): 237/101         -0.129*   -0.137*  
                                 (0.065)   (0.066)  
factor(municip): 240/101         -0.121*   -0.124*  
                                 (0.058)   (0.058)  
factor(municip): 250/101         -0.085    -0.083   
                                 (0.045)   (0.044)  
factor(municip): 251/101          0.067     0.067   
                                 (0.101)   (0.100)  
factor(municip): 253/101         -0.035    -0.035   
                                 (0.033)   (0.032)  
factor(municip): 255/101          0.045     0.037   
                                 (0.058)   (0.058)  
factor(municip): 257/101          0.027     0.025   
                                 (0.083)   (0.083)  
factor(municip): 259/101         -0.068*   -0.064*  
                                 (0.033)   (0.030)  
factor(municip): 260/101         -0.070    -0.054   
                                 (0.054)   (0.052)  
factor(municip): 261/101         -0.075    -0.081   
                                 (0.141)   (0.141)  
factor(municip): 263/101          0.058     0.050   
                                 (0.083)   (0.083)  
factor(municip): 265/101         -0.039    -0.040   
                                 (0.029)   (0.028)  
factor(municip): 267/101          0.050     0.050   
                                 (0.100)   (0.100)  
factor(municip): 269/101         -0.052    -0.054   
                                 (0.048)   (0.047)  
factor(municip): 270/101         -0.038    -0.034   
                                 (0.047)   (0.046)  
factor(municip): 271/101         -0.047    -0.044   
                                 (0.083)   (0.082)  
factor(municip): 301/101         -0.024    -0.009   
                                 (0.116)   (0.114)  
factor(municip): 303/101          0.010     0.032   
                                 (0.116)   (0.115)  
factor(municip): 305/101         -0.120    -0.110   
                                 (0.115)   (0.114)  
factor(municip): 306/101         -0.051    -0.035   
                                 (0.063)   (0.061)  
factor(municip): 307/101         -0.171    -0.162   
                                 (0.101)   (0.099)  
factor(municip): 309/101         -0.087    -0.077   
                                 (0.116)   (0.115)  
factor(municip): 313/101          0.004     0.009   
                                 (0.083)   (0.082)  
factor(municip): 315/101         -0.001    -0.001   
                                 (0.048)   (0.047)  
factor(municip): 316/101         -0.015    -0.007   
                                 (0.045)   (0.043)  
factor(municip): 317/101         -0.087    -0.078   
                                 (0.083)   (0.081)  
factor(municip): 319/101         -0.118    -0.113   
                                 (0.115)   (0.115)  
factor(municip): 320/101         -0.129*   -0.116   
                                 (0.065)   (0.064)  
factor(municip): 321/101         -0.056    -0.057   
                                 (0.090)   (0.090)  
factor(municip): 323/101          0.009     0.012   
                                 (0.052)   (0.050)  
factor(municip): 325/101          0.048     0.072   
                                 (0.053)   (0.051)  
factor(municip): 326/101         -0.089    -0.078   
                                 (0.052)   (0.050)  
factor(municip): 327/101         -0.166*   -0.155*  
                                 (0.077)   (0.076)  
factor(municip): 329/101         -0.072*   -0.070*  
                                 (0.035)   (0.033)  
factor(municip): 330/101         -0.007     0.010   
                                 (0.042)   (0.040)  
factor(municip): 331/101         -0.046    -0.026   
                                 (0.065)   (0.064)  
factor(municip): 333/101         -0.060    -0.048   
                                 (0.042)   (0.040)  
factor(municip): 335/101         -0.061    -0.058   
                                 (0.065)   (0.064)  
factor(municip): 336/101         -0.008     0.000   
                                 (0.058)   (0.056)  
factor(municip): 337/101          0.044     0.048   
                                 (0.198)   (0.197)  
factor(municip): 339/101         -0.151    -0.144   
                                 (0.100)   (0.100)  
factor(municip): 340/101         -0.140*   -0.132*  
                                 (0.055)   (0.053)  
factor(municip): 341/101         -0.053    -0.043   
                                 (0.100)   (0.099)  
factor(municip): 343/101         -0.093    -0.082   
                                 (0.141)   (0.140)  
factor(municip): 345/101         -0.228*   -0.219*  
                                 (0.090)   (0.089)  
factor(municip): 350/101         -0.064    -0.064   
                                 (0.056)   (0.056)  
factor(municip): 351/101          0.003     0.008   
                                 (0.063)   (0.061)  
factor(municip): 353/101          0.053     0.054   
                                 (0.198)   (0.197)  
factor(municip): 355/101         -0.314    -0.287   
                                 (0.198)   (0.198)  
factor(municip): 357/101         -0.095    -0.095   
                                 (0.077)   (0.076)  
factor(municip): 360/101         -0.052    -0.024   
                                 (0.054)   (0.052)  
factor(municip): 361/101          0.037     0.044   
                                 (0.198)   (0.197)  
factor(municip): 363/101          0.103     0.118   
                                 (0.101)   (0.099)  
factor(municip): 365/101         -0.183    -0.171   
                                 (0.116)   (0.115)  
factor(municip): 367/101         -0.190**  -0.168** 
                                 (0.061)   (0.059)  
factor(municip): 369/101         -0.195**  -0.179** 
                                 (0.062)   (0.061)  
factor(municip): 370/101         -0.027    -0.015   
                                 (0.037)   (0.035)  
factor(municip): 371/101          0.003     0.033   
                                 (0.198)   (0.198)  
factor(municip): 373/101         -0.027    -0.023   
                                 (0.038)   (0.037)  
factor(municip): 375/101          0.037     0.049   
                                 (0.090)   (0.089)  
factor(municip): 376/101         -0.035    -0.024   
                                 (0.043)   (0.040)  
factor(municip): 377/101         -0.025    -0.017   
                                 (0.101)   (0.100)  
factor(municip): 379/101         -0.018     0.008   
                                 (0.116)   (0.115)  
factor(municip): 383/101         -0.121    -0.095   
                                 (0.141)   (0.140)  
factor(municip): 385/101          0.027     0.030   
                                 (0.141)   (0.140)  
factor(municip): 387/101         -0.020    -0.005   
                                 (0.090)   (0.089)  
factor(municip): 389/101         -0.216**  -0.203*  
                                 (0.083)   (0.081)  
factor(municip): 390/101          0.024     0.038   
                                 (0.063)   (0.061)  
factor(municip): 391/101         -0.420*   -0.411*  
                                 (0.198)   (0.197)  
factor(municip): 393/101          0.039     0.042   
                                 (0.198)   (0.197)  
factor(municip): 395/101         -0.065    -0.046   
                                 (0.101)   (0.099)  
factor(municip): 397/101         -0.098*   -0.088   
                                 (0.047)   (0.045)  
factor(municip): 400/101         -0.053    -0.023   
                                 (0.043)   (0.042)  
factor(municip): 405/101         -0.048    -0.016   
                                 (0.198)   (0.198)  
factor(municip): 407/101         -0.072    -0.045   
                                 (0.100)   (0.100)  
factor(municip): 410/101         -0.039    -0.040   
                                 (0.058)   (0.057)  
factor(municip): 420/101          0.020     0.031   
                                 (0.052)   (0.050)  
factor(municip): 421/101          0.021     0.041   
                                 (0.083)   (0.081)  
factor(municip): 423/101         -0.019    -0.003   
                                 (0.100)   (0.099)  
factor(municip): 425/101         -0.048    -0.044   
                                 (0.198)   (0.197)  
factor(municip): 427/101         -0.273*** -0.251** 
                                 (0.083)   (0.081)  
factor(municip): 429/101         -0.297**  -0.290** 
                                 (0.100)   (0.100)  
factor(municip): 430/101         -0.126*   -0.118*  
                                 (0.049)   (0.047)  
factor(municip): 431/101         -0.123*   -0.111   
                                 (0.063)   (0.061)  
factor(municip): 433/101          0.047     0.054   
                                 (0.100)   (0.100)  
factor(municip): 435/101          0.099     0.109   
                                 (0.116)   (0.115)  
factor(municip): 437/101          0.043     0.059   
                                 (0.083)   (0.082)  
factor(municip): 439/101         -0.169    -0.168   
                                 (0.116)   (0.115)  
factor(municip): 440/101         -0.062    -0.055   
                                 (0.060)   (0.058)  
factor(municip): 441/101         -0.112    -0.108   
                                 (0.100)   (0.100)  
factor(municip): 443/101         -0.199    -0.181   
                                 (0.198)   (0.198)  
factor(municip): 445/101          0.035     0.039   
                                 (0.054)   (0.053)  
factor(municip): 447/101         -0.135    -0.121   
                                 (0.077)   (0.076)  
factor(municip): 449/101          0.038     0.055   
                                 (0.053)   (0.051)  
factor(municip): 450/101         -0.055    -0.048   
                                 (0.047)   (0.045)  
factor(municip): 451/101         -0.023    -0.024   
                                 (0.100)   (0.100)  
factor(municip): 461/101         -0.025    -0.013   
                                 (0.021)   (0.017)  
factor(municip): 471/101         -0.055    -0.041   
                                 (0.090)   (0.089)  
factor(municip): 473/101         -0.043    -0.037   
                                 (0.065)   (0.064)  
factor(municip): 475/101         -0.213*   -0.171   
                                 (0.101)   (0.101)  
factor(municip): 477/101         -0.053    -0.039   
                                 (0.083)   (0.081)  
factor(municip): 479/101         -0.017    -0.003   
                                 (0.031)   (0.028)  
factor(municip): 480/101         -0.042    -0.030   
                                 (0.072)   (0.071)  
factor(municip): 481/101         -0.083    -0.064   
                                 (0.198)   (0.197)  
factor(municip): 482/101         -0.033    -0.021   
                                 (0.100)   (0.099)  
factor(municip): 483/101         -0.036    -0.036   
                                 (0.101)   (0.100)  
factor(municip): 485/101         -0.030    -0.027   
                                 (0.100)   (0.100)  
factor(municip): 487/101         -0.125    -0.108   
                                 (0.198)   (0.199)  
factor(municip): 489/101         -0.027    -0.023   
                                 (0.090)   (0.089)  
factor(municip): 491/101          0.016     0.027   
                                 (0.100)   (0.099)  
factor(municip): 492/101         -0.106    -0.100   
                                 (0.072)   (0.071)  
factor(municip): 495/101          0.002     0.009   
                                 (0.115)   (0.114)  
factor(municip): 497/101         -0.042    -0.040   
                                 (0.083)   (0.082)  
factor(municip): 499/101          0.291     0.291   
                                 (0.198)   (0.197)  
factor(municip): 501/101         -0.023    -0.027   
                                 (0.116)   (0.115)  
factor(municip): 503/101         -0.075    -0.068   
                                 (0.077)   (0.076)  
factor(municip): 505/101         -0.040    -0.035   
                                 (0.198)   (0.197)  
factor(municip): 507/101         -0.140    -0.138   
                                 (0.101)   (0.100)  
factor(municip): 509/101         -0.065    -0.065   
                                 (0.115)   (0.115)  
factor(municip): 510/101          0.006     0.010   
                                 (0.041)   (0.039)  
factor(municip): 511/101         -0.148    -0.141   
                                 (0.101)   (0.100)  
factor(municip): 513/101         -0.040    -0.037   
                                 (0.101)   (0.100)  
factor(municip): 515/101          0.010     0.024   
                                 (0.047)   (0.045)  
factor(municip): 517/101          0.001     0.017   
                                 (0.141)   (0.140)  
factor(municip): 519/101         -0.048    -0.045   
                                 (0.083)   (0.082)  
factor(municip): 521/101          0.040     0.046   
                                 (0.115)   (0.115)  
factor(municip): 523/101         -0.071    -0.077   
                                 (0.072)   (0.072)  
factor(municip): 525/101         -0.163    -0.157   
                                 (0.141)   (0.140)  
factor(municip): 527/101         -0.128    -0.132   
                                 (0.115)   (0.115)  
factor(municip): 529/101         -0.034    -0.031   
                                 (0.065)   (0.064)  
factor(municip): 530/101         -0.001    -0.002   
                                 (0.056)   (0.056)  
factor(municip): 531/101         -0.018    -0.009   
                                 (0.115)   (0.114)  
factor(municip): 533/101          0.005     0.002   
                                 (0.115)   (0.115)  
factor(municip): 535/101         -0.000    -0.004   
                                 (0.115)   (0.115)  
factor(municip): 537/101         -0.036    -0.026   
                                 (0.048)   (0.046)  
factor(municip): 539/101         -0.047    -0.038   
                                 (0.083)   (0.082)  
factor(municip): 540/101          0.019     0.022   
                                 (0.037)   (0.035)  
factor(municip): 541/101         -0.005     0.007   
                                 (0.072)   (0.071)  
factor(municip): 543/101         -0.058    -0.053   
                                 (0.056)   (0.055)  
factor(municip): 545/101         -0.041    -0.027   
                                 (0.053)   (0.051)  
factor(municip): 550/101          0.010     0.016   
                                 (0.055)   (0.053)  
factor(municip): 551/101         -0.213*   -0.218*  
                                 (0.100)   (0.100)  
factor(municip): 553/101         -0.086    -0.088   
                                 (0.141)   (0.140)  
factor(municip): 555/101          0.145     0.141   
                                 (0.115)   (0.115)  
factor(municip): 557/101         -0.146    -0.149   
                                 (0.083)   (0.083)  
factor(municip): 559/101          0.144     0.144   
                                 (0.115)   (0.115)  
factor(municip): 561/101         -0.022    -0.017   
                                 (0.025)   (0.023)  
factor(municip): 563/101         -0.060    -0.052   
                                 (0.058)   (0.056)  
factor(municip): 565/101         -0.086    -0.087   
                                 (0.072)   (0.072)  
factor(municip): 567/101          0.045     0.043   
                                 (0.116)   (0.115)  
factor(municip): 569/101          0.052     0.053   
                                 (0.090)   (0.090)  
factor(municip): 571/101         -0.112    -0.109   
                                 (0.083)   (0.082)  
factor(municip): 573/101         -0.077    -0.078   
                                 (0.042)   (0.042)  
factor(municip): 575/101         -0.028    -0.026   
                                 (0.036)   (0.035)  
factor(municip): 577/101         -0.049    -0.053   
                                 (0.077)   (0.077)  
factor(municip): 580/101         -0.027    -0.018   
                                 (0.043)   (0.040)  
factor(municip): 601/101         -0.234*   -0.235*  
                                 (0.101)   (0.100)  
factor(municip): 603/101          0.042     0.042   
                                 (0.101)   (0.100)  
factor(municip): 605/101         -0.024    -0.028   
                                 (0.060)   (0.060)  
factor(municip): 607/101          0.001     0.009   
                                 (0.029)   (0.026)  
factor(municip): 609/101         -0.102    -0.101   
                                 (0.115)   (0.115)  
factor(municip): 611/101          0.030     0.030   
                                 (0.090)   (0.090)  
factor(municip): 613/101         -0.053    -0.058   
                                 (0.060)   (0.060)  
factor(municip): 615/101         -0.031    -0.023   
                                 (0.027)   (0.024)  
factor(municip): 617/101          0.078     0.073   
                                 (0.116)   (0.115)  
factor(municip): 619/101         -0.088    -0.089   
                                 (0.072)   (0.072)  
factor(municip): 621/101         -0.059*   -0.057*  
                                 (0.027)   (0.025)  
factor(municip): 623/101         -0.181*   -0.183*  
                                 (0.090)   (0.090)  
factor(municip): 625/101         -0.060    -0.061   
                                 (0.198)   (0.197)  
factor(municip): 627/101         -0.041    -0.044   
                                 (0.101)   (0.100)  
factor(municip): 629/101         -0.010    -0.008   
                                 (0.090)   (0.090)  
factor(municip): 630/101         -0.071    -0.071   
                                 (0.039)   (0.037)  
factor(municip): 631/101         -0.056    -0.048   
                                 (0.038)   (0.036)  
factor(municip): 651/101          0.123     0.121   
                                 (0.101)   (0.100)  
factor(municip): 653/101         -0.053    -0.049   
                                 (0.077)   (0.076)  
factor(municip): 655/101         -0.129    -0.130   
                                 (0.198)   (0.198)  
factor(municip): 657/101         -0.022    -0.019   
                                 (0.028)   (0.026)  
factor(municip): 659/101         -0.168    -0.170   
                                 (0.198)   (0.197)  
factor(municip): 661/101         -0.072*   -0.071*  
                                 (0.032)   (0.031)  
factor(municip): 663/101         -0.062    -0.055   
                                 (0.062)   (0.061)  
factor(municip): 665/101         -0.036    -0.034   
                                 (0.043)   (0.042)  
factor(municip): 667/101          0.080     0.076   
                                 (0.115)   (0.115)  
factor(municip): 669/101         -0.071    -0.074   
                                 (0.115)   (0.115)  
factor(municip): 671/101         -0.028    -0.021   
                                 (0.046)   (0.044)  
factor(municip): 673/101         -0.249*   -0.239*  
                                 (0.115)   (0.114)  
factor(municip): 675/101          0.131     0.135   
                                 (0.141)   (0.140)  
factor(municip): 677/101          0.032     0.027   
                                 (0.141)   (0.141)  
factor(municip): 679/101          0.030     0.034   
                                 (0.090)   (0.090)  
factor(municip): 681/101         -0.035    -0.041   
                                 (0.101)   (0.100)  
factor(municip): 683/101          0.087     0.082   
                                 (0.116)   (0.116)  
factor(municip): 685/101          0.140     0.138   
                                 (0.198)   (0.198)  
factor(municip): 701/101          0.007     0.022   
                                 (0.077)   (0.076)  
factor(municip): 703/101          0.080     0.080   
                                 (0.069)   (0.068)  
factor(municip): 706/101          0.018     0.026   
                                 (0.043)   (0.041)  
factor(municip): 707/101         -0.001     0.015   
                                 (0.045)   (0.043)  
factor(municip): 709/101         -0.078    -0.077   
                                 (0.083)   (0.082)  
factor(municip): 710/101         -0.023    -0.023   
                                 (0.040)   (0.039)  
factor(municip): 711/101          0.009     0.010   
                                 (0.101)   (0.100)  
factor(municip): 713/101         -0.155*   -0.156*  
                                 (0.072)   (0.072)  
factor(municip): 715/101         -0.118    -0.118   
                                 (0.068)   (0.067)  
factor(municip): 717/101          0.026     0.033   
                                 (0.083)   (0.082)  
factor(municip): 719/101         -0.057    -0.050   
                                 (0.141)   (0.140)  
factor(municip): 721/101         -0.013     0.001   
                                 (0.090)   (0.089)  
factor(municip): 723/101          0.005     0.017   
                                 (0.069)   (0.067)  
factor(municip): 725/101         -0.080    -0.064   
                                 (0.141)   (0.140)  
factor(municip): 727/101         -0.017    -0.015   
                                 (0.046)   (0.045)  
factor(municip): 729/101         -0.007    -0.000   
                                 (0.083)   (0.082)  
factor(municip): 730/101         -0.058    -0.050   
                                 (0.031)   (0.029)  
factor(municip): 731/101         -0.042    -0.030   
                                 (0.032)   (0.029)  
factor(municip): 733/101         -0.018    -0.016   
                                 (0.083)   (0.082)  
factor(municip): 735/101         -0.052    -0.037   
                                 (0.116)   (0.115)  
factor(municip): 737/101          0.007     0.008   
                                 (0.060)   (0.059)  
factor(municip): 739/101         -0.020    -0.014   
                                 (0.077)   (0.076)  
factor(municip): 740/101          0.021     0.024   
                                 (0.033)   (0.031)  
factor(municip): 741/101         -0.137    -0.119   
                                 (0.198)   (0.197)  
factor(municip): 743/101         -0.036    -0.031   
                                 (0.033)   (0.031)  
factor(municip): 745/101         -0.119*   -0.119*  
                                 (0.060)   (0.059)  
factor(municip): 746/101         -0.011    -0.012   
                                 (0.041)   (0.040)  
factor(municip): 747/101          0.012     0.020   
                                 (0.083)   (0.082)  
factor(municip): 749/101         -0.093    -0.091   
                                 (0.100)   (0.100)  
factor(municip): 751/101         -0.025    -0.018   
                                 (0.018)   (0.015)  
factor(municip): 756/101         -0.097    -0.088   
                                 (0.053)   (0.051)  
factor(municip): 760/101         -0.010    -0.006   
                                 (0.043)   (0.041)  
factor(municip): 761/101          0.013     0.003   
                                 (0.141)   (0.141)  
factor(municip): 763/101         -0.092    -0.087   
                                 (0.090)   (0.089)  
factor(municip): 765/101         -0.211    -0.208   
                                 (0.198)   (0.198)  
factor(municip): 766/101          0.032     0.031   
                                 (0.048)   (0.048)  
factor(municip): 767/101          0.124     0.122   
                                 (0.090)   (0.090)  
factor(municip): 769/101         -0.108    -0.104   
                                 (0.198)   (0.197)  
factor(municip): 771/101          0.005     0.006   
                                 (0.077)   (0.076)  
factor(municip): 773/101         -0.094*   -0.084   
                                 (0.047)   (0.045)  
factor(municip): 775/101          0.155     0.155   
                                 (0.198)   (0.197)  
factor(municip): 777/101         -0.198*   -0.189*  
                                 (0.090)   (0.089)  
factor(municip): 779/101         -0.067    -0.061   
                                 (0.035)   (0.033)  
factor(municip): 781/101          0.054     0.056   
                                 (0.141)   (0.140)  
factor(municip): 783/101         -0.031    -0.024   
                                 (0.115)   (0.115)  
factor(municip): 785/101          0.010     0.007   
                                 (0.072)   (0.072)  
factor(municip): 787/101         -0.048    -0.046   
                                 (0.036)   (0.035)  
factor(municip): 789/101          0.070     0.076   
                                 (0.198)   (0.198)  
factor(municip): 791/101         -0.009    -0.007   
                                 (0.030)   (0.029)  
factor(municip): 793/101         -0.062    -0.054   
                                 (0.090)   (0.089)  
factor(municip): 801/101         -0.070    -0.062   
                                 (0.141)   (0.140)  
factor(municip): 803/101         -0.013     0.015   
                                 (0.198)   (0.197)  
factor(municip): 805/101         -0.031    -0.015   
                                 (0.049)   (0.047)  
factor(municip): 807/101         -0.076    -0.063   
                                 (0.065)   (0.064)  
factor(municip): 809/101          0.180     0.179   
                                 (0.198)   (0.197)  
factor(municip): 810/101         -0.159*** -0.150***
                                 (0.046)   (0.044)  
factor(municip): 811/101         -0.125    -0.115   
                                 (0.065)   (0.064)  
factor(municip): 813/101         -0.028    -0.004   
                                 (0.031)   (0.029)  
factor(municip): 815/101         -0.088    -0.070   
                                 (0.100)   (0.099)  
factor(municip): 817/101         -0.043    -0.012   
                                 (0.115)   (0.114)  
factor(municip): 819/101         -0.020     0.005   
                                 (0.069)   (0.068)  
factor(municip): 820/101         -0.009    -0.000   
                                 (0.062)   (0.061)  
factor(municip): 821/101         -0.005     0.009   
                                 (0.051)   (0.050)  
factor(municip): 823/101         -0.055    -0.045   
                                 (0.115)   (0.115)  
factor(municip): 827/101         -0.334*** -0.322** 
                                 (0.100)   (0.099)  
factor(municip): 829/101         -0.135    -0.113   
                                 (0.100)   (0.100)  
factor(municip): 831/101         -0.027    -0.022   
                                 (0.083)   (0.081)  
factor(municip): 833/101         -0.133    -0.128   
                                 (0.116)   (0.115)  
factor(municip): 835/101         -0.077    -0.045   
                                 (0.115)   (0.116)  
factor(municip): 837/101          0.010     0.023   
                                 (0.072)   (0.071)  
factor(municip): 839/101          0.115     0.137   
                                 (0.090)   (0.089)  
factor(municip): 840/101         -0.054    -0.051   
                                 (0.056)   (0.055)  
factor(municip): 841/101         -0.014     0.016   
                                 (0.062)   (0.062)  
factor(municip): 843/101         -0.053    -0.050   
                                 (0.141)   (0.140)  
factor(municip): 845/101         -0.093    -0.095   
                                 (0.101)   (0.100)  
factor(municip): 846/101         -0.038    -0.030   
                                 (0.052)   (0.050)  
factor(municip): 847/101          0.027     0.042   
                                 (0.058)   (0.056)  
factor(municip): 849/101         -0.059    -0.049   
                                 (0.039)   (0.036)  
factor(municip): 851/101         -0.027    -0.013   
                                 (0.021)   (0.017)  
factor(municip): 860/101         -0.017    -0.002   
                                 (0.036)   (0.033)  
factor(municip): 861/101         -0.031    -0.026   
                                 (0.062)   (0.061)  
factor(essround): 2/1             0.014     0.021*  
                                 (0.009)   (0.010)  
factor(essround): 3/1             0.064***  0.059***
                                 (0.009)   (0.009)  
factor(essround): 4/1            -0.078*** -0.094***
                                 (0.010)   (0.015)  
factor(essround): 5/1            -0.152*** -0.148***
                                 (0.009)   (0.010)  
Unemployed                       -0.051*** -0.050***
                                 (0.010)   (0.010)  
Age                               0.000     0.000   
                                 (0.000)   (0.000)  
Years of education                0.002*    0.002*  
                                 (0.001)   (0.001)  
Male                              0.046***  0.046***
                                 (0.005)   (0.005)  
Right--wing                       0.165***  0.166***
                                 (0.013)   (0.013)  
Political interest               -0.008*   -0.008*  
                                 (0.004)   (0.004)  
Household Unemployment           -0.021*   -0.021*  
                                 (0.011)   (0.011)  
Income (in 100k DKR)              0.004*    0.004*  
                                 (0.002)   (0.002)  
Years lived at location           0.000     0.000   
                                 (0.000)   (0.000)  
Population (80m)                 -0.000    -0.000** 
                                 (0.000)   (0.000)  
Local Unemployment Rate (80m)    -0.112*   -0.099*  
                                 (0.046)   (0.044)  
Newspaper exposure                0.087     0.078   
                                 (0.059)   (0.043)  
popdens2500                      -0.000             
                                 (0.000)            
Local Unemployment Rate (2500m)   0.245             
                                 (0.142)            
ueb80 x nwsptot                   0.230     0.197   
                                 (0.172)   (0.165)  
nwsptot x ueb2500                -0.959*            
                                 (0.478)            
munemploy                                  -0.066   
                                           (0.455)  
nwsptot x munemploy                        -2.245** 
                                           (0.800)  
----------------------------------------------------
R-squared                            0.219     0.220
adj. R-squared                       0.174     0.174
sigma                                0.197     0.197
F                                    4.790     4.821
p                                    0.000     0.000
Log-likelihood                    1415.067  1417.234
Deviance                           218.762   218.604
AIC                              -2162.135 -2168.468
BIC                                 75.057    62.026
N                                 5992      5992    
====================================================

> #F-tests for interaction terms in FE models
> fit20<-lm(stfeco~factor(municip)+factor(essround)+unemplo5+ageR+eduresu+maleR+lrscale+polintr
+ 	+pincome2000_1k+dayslived_y+popdens80+ueb80+munemploy+nwsptot
+ 	,data=dat[rownames(dat)%in%rownames(fit2FE$model),])

> fit0<-lm(stfeco~factor(municip)+factor(essround)+unemplo5+ageR+eduresu+maleR+lrscale+polintr
+ 	+pincome2000_1k+dayslived_y+popdens80+ueb80+nwsptot+popdens2500+ueb2500
+ 	,data=dat[rownames(dat)%in%rownames(fit2FE$model),])

> anova(fit0,fit1FE)
Analysis of Variance Table

Model 1: stfeco ~ factor(municip) + factor(essround) + unemplo5 + ageR + 
    eduresu + maleR + lrscale + polintr + pincome2000_1k + dayslived_y + 
    popdens80 + ueb80 + nwsptot + popdens2500 + ueb2500
Model 2: stfeco ~ factor(municip) + factor(essround) + unemplo4 + ageR + 
    eduresu + maleR + lrscale + polintr + oneunemplo4 + pincome2000_1k + 
    dayslived_y + popdens80 + ueb80 * nwsptot + popdens2500 + 
    ueb2500 * nwsptot
  Res.Df    RSS Df Sum of Sq      F  Pr(>F)  
1   5662 219.18                              
2   5659 218.76  3   0.42111 3.6312 0.01237 *
---
Signif. codes:  0 �***� 0.001 �**� 0.01 �*� 0.05 �.� 0.1 � � 1

> anova(fit20,fit2FE)
Analysis of Variance Table

Model 1: stfeco ~ factor(municip) + factor(essround) + unemplo5 + ageR + 
    eduresu + maleR + lrscale + polintr + pincome2000_1k + dayslived_y + 
    popdens80 + ueb80 + munemploy + nwsptot
Model 2: stfeco ~ factor(municip) + factor(essround) + unemplo4 + ageR + 
    eduresu + maleR + lrscale + polintr + oneunemplo4 + pincome2000_1k + 
    dayslived_y + popdens80 + ueb80 * nwsptot + munemploy * nwsptot
  Res.Df    RSS Df Sum of Sq     F   Pr(>F)   
1   5663 219.17                               
2   5660 218.60  3   0.56822 4.904 0.002097 **
---
Signif. codes:  0 �***� 0.001 �**� 0.01 �*� 0.05 �.� 0.1 � � 1

> ######################################
> # Code for Figure 3			#
> 
> summary(dat$munemploy)
   Min. 1st Qu.  Median    Mean 3rd Qu.    Max.    NA's 
0.00800 0.03600 0.05000 0.05031 0.06400 0.15500     150 

> news<-seq(0,1,length.out=10)

> muni<-seq(min(dat$munemploy,na.rm=TRUE),max(dat$munemploy,na.rm=TRUE),length.out=10)

> ueb25<-seq(min(dat$ueb2500,na.rm=TRUE),max(dat$ueb2500,na.rm=TRUE),length.out=10)

> fit1<-fit1FE #figure was originally without municipality FEs and this adjusts the code to models /w FEs

> fit2<-fit2FE

> #coefs
> b1<-fit1$coef["ueb2500"]

> b3<-fit1$coef["nwsptot:ueb2500"]

> #SEs
> v1<-diag(vcov(fit1))["ueb2500"]

> v3<-diag(vcov(fit1))["nwsptot:ueb2500"]

> cov3<-vcov(fit1)["nwsptot:ueb2500","ueb2500"]

> se<-rep(NA,10)

> cf<-rep(NA,10)

> for(i in 1:10){
+ 	se[i]<-sqrt(v1+news[i]^2*v3+2*news[i]*cov3)
+ 	cf[i]<-b1+b3*news[i]
+ 	}

> se
 [1] 0.1419185 0.1134805 0.1061053 0.1236221 0.1579558 0.2006520 0.2474186
 [8] 0.2963347 0.3464912 0.3974186

> cf
 [1]  0.24494688  0.13842734  0.03190780 -0.07461174 -0.18113128 -0.28765082
 [7] -0.39417036 -0.50068991 -0.60720945 -0.71372899

> w<-matrix(NA,nr=length(ueb.names),nc=2)#weighting wrt to variation (1st to 99th percentile)

> for (i in 1:length(ueb.names))w[i,]<-quantile(dat[dat$essround!=9&dat$popdens80>9,ueb.names[i]],c(.01,.99),na.rm=T)

> w.out<-apply(w,1,diff)

> se<-se*w.out[14]	

> cf<-cf*w.out[14]

> #coefs
> b1<-fit2$coef["munemploy"]

> b3<-fit2$coef["nwsptot:munemploy"]

> #SEs
> v1<-diag(vcov(fit2))["munemploy"]

> v3<-diag(vcov(fit2))["nwsptot:munemploy"]

> cov3<-vcov(fit2)["nwsptot:munemploy","munemploy"]

> se1<-rep(NA,10)

> cf1<-rep(NA,10)

> for(i in 1:10){
+ 	se1[i]<-sqrt(v1+news[i]^2*v3+2*news[i]*cov3)
+ 	cf1[i]<-b1+b3*news[i]
+ 	}

> se1
 [1] 0.4545382 0.4316465 0.4264316 0.4395230 0.4693916 0.5131158 0.5675021
 [8] 0.6297943 0.6978785 0.7702204

> cf1
 [1] -0.06616358 -0.31556775 -0.56497193 -0.81437611 -1.06378028 -1.31318446
 [7] -1.56258863 -1.81199281 -2.06139699 -2.31080116

> w.mun<-quantile(dat[dat$essround!=9&dat$popdens80>9,"munemploy"],c(0.01,.99),na.rm=T)

> w.mun<-w.mun[2]-w.mun[1]

> se1<-se1*w.mun

> cf1<-cf1*w.mun

> #coefs
> b1<-fit2$coef["ueb80"]

> b3<-fit2$coef["ueb80:nwsptot"]

> #SEs
> v1<-diag(vcov(fit2))["ueb80"]

> v3<-diag(vcov(fit2))["ueb80:nwsptot"]

> cov3<-vcov(fit2)["ueb80:nwsptot","ueb80"]

> se2<-rep(NA,10)

> cf2<-rep(NA,10)

> for(i in 1:10){
+ 	se2[i]<-sqrt(v1+news[i]^2*v3+2*news[i]*cov3)
+ 	cf2[i]<-b1+b3*news[i]
+ 	}

> se2
 [1] 0.04373112 0.03333803 0.03132204 0.03900447 0.05227102 0.06792498
 [7] 0.08465210 0.10192540 0.11950828 0.13728183

> cf2
 [1] -0.09908883 -0.07715785 -0.05522687 -0.03329589 -0.01136491  0.01056607
 [7]  0.03249705  0.05442803  0.07635901  0.09828999

> se2<-w.out[1]*se2

> cf2<-w.out[1]*cf2

> #######PLOTTING####
> #pdf('nwsp-int.pdf',width=8,height=3)
> setEPS()

> postscript('nwsp-int.eps',width=8,height=3)

> par(mfrow=c(1,3),omd=c(.05,1,0,1))

> #par(mfrow=c(1,1))
> plot(news,cf2,type="l",ylim=c(-.2,.1),cex.axis=.8,,main="",
+ 	ylab="",xaxt="n",yaxt="n",
+ 	xlab="")

> axis(2,at=seq(-.2,.1,.05),las=2,cex.axis=1)

> axis(1,cex.axis=1)

> mtext(2,text="Marginal effect of unemployment rate",line=5)

> mtext(2,text="(1st to 99th percentile)",line=3.75,cex=.8)

> mtext(3,text="A",font=2,line=2,cex=1.5)

> mtext(3,text="Unemployment in 80m radius",font=1,line=0.5,cex=.8)

> abline(h=0,lty=2)

> #X.vec<-c(news,1,rev(news),0)
> #Y.vec<-c(cf-1.97*se,cf[10]+1.97*se[10],rev(cf+1.97*se),cf[1]+1.97*se[1])
> #polygon(X.vec,Y.vec,col="grey90",border=NA);lines(news,cf,col="black")
> 
> lines(news,cf2+1.97*se2)

> lines(news,cf2-1.97*se2)

> par(new=TRUE)

> plot.new()

> hist(dat$nwsptot[dat$popdens80>9&dat$essround!=9],xpd=TRUE,col="grey95",breaks=na.omit(unique(dat$nwsptot)),ylim=c(0,15),
+ 	freq=FALSE,axes=FALSE,xlab="",ylab="",main="")

> #dev.new(width=4,height=4.7)
> #par(mfrow=c(1,1))
> plot(news,cf,type="l",ylim=c(-.2,.1),cex.axis=.8,,main="",
+ 	ylab="",xaxt="n",yaxt="n",
+ 	xlab="")

> axis(2,at=seq(-.2,.1,.05),las=2,cex.axis=1)

> axis(1,cex.axis=1)

> mtext(3,text="B",font=2,line=2,cex=1.5)

> mtext(3,text="Unemployment in 2500m radius",font=1,line=0.5,cex=.8)

> mtext(1,text="Newspaper exposure",line=3)

> abline(h=0,lty=2)

> #X.vec<-c(news,1,rev(news),0)
> #Y.vec<-c(cf-1.97*se,cf[10]+1.97*se[10],rev(cf+1.97*se),cf[1]+1.97*se[1])
> #polygon(X.vec,Y.vec,col="grey90",border=NA);lines(news,cf,col="black")
> 
> lines(news,cf+1.97*se)

> lines(news,cf-1.97*se)

> par(new=TRUE)

> plot.new()

> hist(dat$nwsptot[dat$popdens80>9&dat$essround!=9],xpd=TRUE,col="grey95",breaks=na.omit(unique(dat$nwsptot)),ylim=c(0,15),
+ 	freq=FALSE,axes=FALSE,xlab="",ylab="",main="")

> #dev.new(width=4,height=4.7)
> #par(fig=c(0,1,0,1))
> plot(news,cf1,type="l",ylim=c(-.2,.1),cex.axis=.8,main="",
+ 	ylab="",xaxt="n",yaxt="n",
+ 	xlab="")

> axis(2,at=seq(-.2,.1,.05),las=2,cex.axis=1)

> axis(1,cex.axis=1)

> mtext(3,text="C",font=2,line=2,cex=1.5)

> mtext(3,text="Unemployment in municipality",font=1,line=0.5,cex=.8)

> lines(news,cf1+1.97*se1)

> lines(news,cf1-1.97*se1)

> abline(h=0,lty=2)

> par(new=TRUE)

> plot.new()

> hist(dat$nwsptot[dat$popdens80>9&dat$essround!=9],xpd=TRUE,col="grey95",breaks=na.omit(unique(dat$nwsptot)),ylim=c(0,15),
+ 	freq=FALSE,axes=FALSE,xlab="",ylab="",main="")

> dev.off()
null device 
          1 

> #END OF SCRIPT
